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Abstract 

The quark (gluon)-hadron duality constitutes a basis for the theoretical 
treatment of a wide range of inclusive processes - from hadronic r decays and 
R e + e - , to semileptonic and nonleptonic decay rates of heavy flavor hadrons. 
A theoretical analysis of these processes is carried out by using the operator 
product expansion (OPE) in the Euclidean domain, with subsequent analytic 
continuation to the Minkowski domain. We formulate the notion of the quark 
(gluon)-hadron duality in quantitative terms, then classify various contribu- 
tions leading to violations of duality. A prominent role in the violations of 
duality seems to belong to the so called exponential terms which, conceptually, 
may represent the (truncated) tail of the power series. A qualitative model, 
relying on an instanton background field, is developed, allowing one to get an 
estimate of the exponential terms. We then discuss a number of applications, 
mostly from heavy quark physics. 



CERN-TH/96-113 
April 1996 



1 Introduction 



Nonperturbative effects have been analyzed in QCD in the framework of the operator 
product expansion (OPE) Jl], || since the inception of QCD. Recently a remarkable 
progress has been achieved along these lines in the heavy quark theory (for reviews 
see Ref. ||). A large number of applications of OPE in heavy quark theory refer 
to quantities of the essentially Minkowski nature, e.g. calculations of the inclusive 
decay widths, spectra and so on. Wilson's operator product expansion per se is 
formulated in the Euclidean domain. The expansion built in the Euclidean domain 
and by necessity truncated, is translated in the language of the observables through 
an analytic continuation. An indispensable element of this procedure, the so called 
quark-hadron duality, is always assumed, most often tacitly. This paper is devoted 
to the discussion of quark-hadron duality and deviations from it. Although the issue 
will be considered primarily in the context of heavy quark theory, the problems we 
will deal with are quite general and are by no means confined to heavy quark theory. 
For instance, determination of a s from the hadronic width of r - a problem of 
paramount importance now under intense scrutiny M - falls into this category. It is 
known || that deviations from duality may be conceptually related to the behavior 
of the operators of high dimension in OPE. Unfortunately, very little is known 
about this behavior in the quantitative sense, beyond the fact that the expansion is 
asymptotic |BJ. Therefore we are forced to approach the problem from the other side 
- engineering a model which allows us to start discussing deviations from duality. 
The model is based on instantons, but by no means is derivable in QCD. Moreover, it 
does not exhaust all mechanisms which might lead, in principle, to deviations from 
duality, focusing, rather, on one specific contribution - the so called exponential 
terms. Nevertheless, it seems to be physically motivated and can serve for qualitative 
analysis at present, and as a guideline for future refinements. 

Indeed, the (fixed size) instanton contribution to correlation functions with large 
momentum transfers can be interpreted as a mechanism in which the large external 
momentum is transmitted through a soft coherent field configuration. Speaking 
graphically, the large external momentum is shared by a very large number of quanta 
so that each quantum is still relatively soft. It is clear that this mechanism is not 
represented in the practical version of OPE |2[], and, thus, gives an idea of how 
strong deviations from duality might be. 

One of the most interesting aspects revealed in this model is the distinct nature 
of exponential contributions absent in practical OPE, both in the Euclidean and the 
Minkowski domains. If in the former the exponential effects die off fast enough, in 
the latter, deviations from duality are suppressed to a lesser extent - the exponential 
fall off is milder, and it is modulated by oscillations. These features seem to be so 
general that most certainly they will survive in future treatments which, hopefully, 
will be significantly closer to fundamental QCD than our present consideration. 

If one accepts this model, at least for orientation, many interesting technical 
problems arise. Instanton contributions in heavy quark theory were previously dis- 
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cussed more than once || [7]. Although the corresponding analysis seemed rather 
straightforward at first sight, it resulted in some apparent paradoxes; for instance, 
the instanton contribution to the spectrum of the inclusive heavy quark decays 
seemingly turned out to be parametrically larger than the very same contribution 
to the total decay rate JF[ ||. The puzzle is readily solvable, however: one observes 
that the problem lies in the separation of the exponentially small terms from the 
"background" of the power terms of OPE; this is a subtle and, generally speaking, 
ambiguous procedure, particularly in the Minkowski domain, and depends on the 
specific quantity under consideration. We will dwell on this issue at length in the 
present paper. 

We begin, however, with a brief formulation of the very notion of duality (Sect. 
2). Quantifying this notion is an important task by itself. In modeling deviations 
from duality, the adoption of the following attitude is made so as to stay on safe 
ground: we will try to develop a model yielding a conservative estimate on the upper 
bound for deviations from duality. In other words, given the prediction for this or 
that quantity based on duality (i.e. spectra, total inclusive widths and so on), one 
establishes the accuracy with which this prediction is expected to be valid. In this 
way one sets the lower limit on the energy release needed to achieve the required 
accuracy. For this limited purpose even a crude model, such as the instanton model 
to be discussed below, may be sufficient, perhaps, after some minor refinements. 

Why is this attitude logical? If we knew in detail some specific mechanism 
omitted in the theoretical calculations - whether associated with the truncation 
of practical OPE, or due to other sources - we could include it in the theoretical 
prediction for the cross sections and say that the actual hadronic cross section is 
dual to this new improved prediction. Thus, paradoxically, the very nature of duality 
implies that deviations from it are always estimated roughly. Analyzing deviations 
from duality at each given stage of development of QCD is equivalent to analyzing 
our ignorance, rather than our knowledge. At the present stage, as was already 
mentioned, our knowledge is, more or less, limited to practical OPE. 

A much more ambitious goal is developing a framework suitable for actual cal- 
culation of extra contributions not seen in practical OPE. Although the instanton 
model is sometimes used for this purpose as well, one should clearly realize that 
quantitatively reliable results are not expected to emerge in this way. This is a 
speculative procedure intended only for qualitative orientation. We will occasion- 
ally resort to it only due to the absence of better ideas. One may hope that a 
universal qualitative picture will be revealed en route, which will be robust enough 
to survive future developments of the issue. 

Although this problem - estimates of deviations from duality - is obviously of 
paramount practical importance, surprisingly little has been said about this subject 
in the literature. Apart from some general remarks presented in Ref. ||, an attempt 
to discuss the issue in a different (exclusive) context was made in @. 

Our paper is organized below as follows: In Sect. 2 we outline the general 
principles behind duality and its violation. Sect. 3 is devoted to general features of 
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the exponential terms believed to be responsible for duality violation. In particular, 
the distinction between their patterns in the Euclidean and Minkowski domains is 
explained here. In Sect. 4 we outline the instanton model we use as a framework to 
generate exponential terms. Sect. 5 illustrates our main points in what is, probably, 
the most transparent example: e + e~ annihilation and the hadronic decays of the 
r lepton. In Sect. 6 we discuss the general features of heavy quark decays in the 
instanton background. In Sect. 7, we begin the business of actual calculation - to 
warm up, we consider a toy model where the spins of all relevant fields are discarded 
to avoid technicalities. Section 8 is devoted to actual heavy quark decays in QCD. 
The exponential contributions are estimated, both in the spectra, and in the inclusive 
decay rates, for the transitions of the heavy quark into a massless one. In section 9 we 
address the applied, but practically important, problem of deviations from duality 
in the semileptonic decays of D and B mesons. There are good phenomenological 
reasons to believe that in D decays these deviations are significant, of order 0.5. 
Adjusting parameters of the model in such a way as to explain these deviations, we 
conclude that deviations from duality in the B decays are expected to be negligibly 
small (in the total semileptonic decay rate and in the similar radiative processes). 
The effect seems to be larger - perhaps even detectable - in the inclusive nonleptonic 
rates. The drawbacks and deficiencies of the model we use for the estimates of the 
exponential terms are summarized in Sect. 10. We present some comments on 
the vast literature treating the processes under discussion in Sect. 11. Section 12 
summarizes our results and outlines problems for future exploration. 

2 Duality and the OPE 

Wilson's OPE is the basis of virtually all calculations of nonperturbative effects in 
analytical QCD. Since the very definition of duality relies heavily on Wilson's OPE, 
we first briefly review its main elements. For the sake of definiteness, we will speak 
of the heavy quark expansion, although one should keep in mind that the procedure 
is quite general; in other processes (e.g. the hadronic r decays) the wording must 
be somewhat changed, but the essence remains intact. 

The original QCD Lagrangian is formulated at very short distances. Starting 
from this Lagrangian, one evolves it down, integrating out all fluctuations with 
frequencies /i < u < Mo where Mo is the original normalization point, and /i will 
be treated, for the time being, as a current parameter. In this way we get the 
Lagrangian which has the form 

C = J2C n (M ;^)OM. (1) 

n 

The coefficient functions, C n , represent the contribution of virtual momenta from 
H to M . The operators, O n , enjoy the full rights of Heisenberg operators with 
respect to all field fluctuations with frequencies less than \i. The sum in Eq. (|TJ) is 
infinite - it runs over all possible Lorentz singlet gauge invariant operators which 
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possess the appropriate quantum numbers. The operators can be ordered according 
to their dimension; moreover, we can use the equations of motion, stemming from 
the original QCD Lagrangian, to get rid of some of the operators in the sum. Those 
operators that are reducible to full derivatives give vanishing contributions to the 
physical (on mass shell) matrix elements, and can thus be discarded as well. 

Speaking abstractly, one is free to take any value of /i in Eq. (|I|); in particular, 
fi = would mean that everything is calculated and we have the full 5* matrix, all 
conceivable amplitudes, at our disposal. Nothing is left to be done. In this case 
Eq. (P is just a sum of all possible amplitudes. This sum then must be written 
in terms of physical hadronic states, of course, not in terms of the quark and gluon 
operators since the latter degrees of freedom simply do not survive scales below some 

Needless to say, present-day QCD does not allow the explicit evolution down to 
fi = 0. Calculating the coefficient functions we have to stop somewhere, at such 
virtualities that the quark and gluon degrees of freedom are still relevant, and the 
coefficient functions C n (Mo, /i) are still explicitly calculable. On the other hand, for 
obvious reasons, it is highly desirable to have fi as low as possible. In the heavy 
quark theory there is an additional requirement that [i must be much less than tuq. 

Let us assume that /i is large enough so that a s (fi)/n <C 1 on the one hand, and 
small enough so that there is no large gap between Aqcd and /x. The possibility 
to make such a choice of /i could not be anticipated a priori and is an extremely 
fortunate feature of QCD. Quarks and gluons with offshellness larger than [i chosen 
in this way will be referred to as hard. 

All observable amplitudes must be /i independent, of course. The \i dependence 
of the coefficient functions C n must conspire with that of the matrix elements of 
the operators O n in such a way as to ensure this /i independence of the physical 
amplitudes. 

What can be said about the calculation of the coefficients C n ? Since /i is suffi- 
ciently large, the main contribution comes from perturbation theory. We just draw 
all relevant Feynman graphs and calculate them, generating an expansion in a s (/i) 

C n = X)aja£(/x) . (2) 
i 

(Sometimes some graphs will contain not only the powers of ce s (n) but also powers 
of a s ln(rriQ/ [i). This happens if the anomalous dimension of the operator O n is 
nonvanishing, or if a part of a contribution to C n comes from characteristic momenta 
of order ttlq and is, thus, expressible in terms of a s (mg), and we rewrite it in terms 
of a s (/i).) 

As a matter of fact expression is not quite accurate theoretically. One should 
not forget that, in doing the loop integrations in C n , we must discard the domain 
of virtual momenta below fi, by definition of C n (fi). Subtracting this domain from 
the perturbative loop integrals, we introduce in C n power corrections of the type 
{n/mQ) n by hand. In principle, one should recognize the existence of such corrections 
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and deal with them. The fact that they are actually present was realized long ago @ . 
Neglecting them, at the theoretical level, results in countless paradoxes which still 
surface from time to time in the literature. If it is possible to choose \i sufficiently 
small, these corrections may be insignificant numerically, and can be omitted. This 
is what is actually done in practice. This is one of the elements of a simplification of 
the Wilson operator product expansion. The simplified version is called the practical 
version of the OPE, or practical OPE ||. 

Even if perturbation theory may dominate the coefficient functions, they still 
also contain nonperturbative terms coming from short distances. Sometimes they 
are referred to as noncondensate nonperturbative terms. An example is provided 
by direct instantons flIU| , with sizes of order uiq . These contributions fall off as 
high powers of Aq CD /mQ (or Aq CD /E where E is a characteristic energy release in 
the process under consideration), and are very poorly controlled theoretically. Since 
the fall off of the noncondensate nonperturbative corrections is extremely steep, 
basically the only thing we can say about them is that there is a critical value of 
rriQ (or E). For lower values of tuq (or E) no reliable theoretical predictions are 
possible at present. For higher values of ttlq one can ignore the noncondensate 
nonperturbative contributions. The noncondensate nonperturbative contributions 
are neglected in the practical OPE. In what follows, we will not touch upon these 
type of effects which are associated with the (small-size) instanton contributions 
to the coefficient functions. There is another, technical, reason why we choose 
not to consider these effects. Since the small-size instantons represent hard field 
fluctuations, all heavy quark expansions carried out in the spirit of HQET (Tl| 
become invalid; the corresponding theory has to be developed anew. In particular, 
the standard HQET decomposition of the heavy quark field in the form Q(x) = 
exp{imQV /i a; At }Q(x) becomes inapplicable, as well as the statement that all heavy 
quark spin effects are suppressed by 1/rriQ, and so on. This circumstance is not fully 
recognized in the literature. Due to these reasons, we instead focus on effects due 
to large-size instantons. This will provide a workable framework for visualizing the 
exponential term. 

At very large tuq (or E), the exponential terms are parametrically smaller (in the 
Euclidean domain) than the power-like non-condensate nonperturbative corrections 
in the coefficient functions. One can argue, however, that this natural hierarchy 
sets in at such large values of momentum transfer where both effects are practically 
unimportant. At intermediate values of the momentum transfers - most interesting 
from the point of view of applications - an inverse hierarchy may take place, where 
the exponential terms are numerically more important. 

Ignoring the nonperturbative contributions in the coefficient functions is not the 
only simplification in the practical OPE. The series of operators appearing in £ 
(the condensate series) is infinite. Practically we truncate it in some finite order, so 
that the sum in the expansion we deal with approximates the exact result, but by 
no means coincides with it. The truncation of the expansion is a key point. The 
condensate expansion is asymptotic |J. Therefore, expanding it to higher orders 
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indefinitely does not mean that the accuracy of the approximation to the exact 
result becomes better. On the contrary, as in any asymptotic series, there exists an 
optimal order. Truncating the series at this order, we get the best accuracy. The 
difference between the exact result and the series truncated at the optimal order is 
exponential. Large-size instantons, treated in an appropriate way will, in a sense, 
represent the high-order tail omitted in the truncated series. 

The essence of the phenomenon - occurrence of the exponential terms - is similar 
to the emergence of the condensates at a previous stage. Indeed, let us consider, first, 
the standard Feynman perturbation theory. At any finite order the perturbative 
contribution is well-defined. At the same time, the coefficients of the a s series grow 
factorially with n, and this means that the a s series must be, somehow, cut off, i.e. 
regularized. The proper way of handling this factorial divergence is by introducing 
the normalization point and the condensate corrections which tempers the factorial 
divergence of the Feynman perturbative series in high orders and, simultaneously, 
bring in terms of order exp (—C/a s {mQ)) where C is some positive constant. Loosely 
speaking, one may say that contributions of this type are related to the high-order 
tails of the a s series. Similarly, the high-order tails of the condensate (power) series 
correspond to the occurrence of the exponential terms. Correspondingly, the OPE, 
even optimally truncated, approximates the exact result up to exponential terms. 

The exponential terms not seen in the practical OPE appear both in Euclidean, 
and Minkowski quantities. Their particular roles and behaviors are quite different, 
however. Technically, the rate of fall off is much faster in the Euclidean domain 
than in the Minkowski domain, as we will see later. Conceptually, the exponential 
terms in the Minkowski domain determine deviations from duality. 

Let us now describe what we mean by duality in somewhat more detail. Assume 
that the effective Lagrangian we work with includes external sources, so that the 
expectation value of this Lagrangian actually yields the complete set of physical 
amplitudes. The physically observable Minkowski quantities (i.e. spectra, total 
hadronic widths and so on) are given by the imaginary parts of certain terms in the 
effective Lagrangian. These terms are calculated as an expansion in the Euclidean 
domain. This is a practical necessity - since our theoretical tools are based on 
the expansions phrased in terms of quarks and gluons, we have to operate in the 
Euclidean domain. We then analytically continue in relevant momentum transfers to 
the Minkowski domain. Of course, if we could find the exact result in the Euclidean 
domain, its analytic continuation to the Minkowski domain would yield the exact 
spectra, etc., - there would be no need in introducing duality at all. In reality, 
the calculation is done using the practical OPE. Both, the perturbative series in 
the coefficient functions and the condensate series are truncated at a certain order. 
We then analytically continue each individual term in the expansion thus obtained, 
term by term, from the deep Euclidean domain to the Minkowski one, and take the 
imaginary part. The corresponding prediction, which can be interpreted in terms 
of quarks and gluons, is declared to be dual to physically measurable quantities in 
terms of hadrons provided that the energy release is large. In this context dual means 



6 



approximately equal. The discrepancy between the exact (hadronic) result and the 
quark- gluon prediction based on the practical OPE is referred to as a deviation from 
(local) duality. 

Reiterating, in defining duality, we first do a straightforward analytical continu- 
ation, term by term. Defining the analytic continuation to the Minkowski domain 
in this way, it is not difficult to see, diagrammatically, that those lines which were 
far off shell in the Euclidean calculation remain hard in the sense that now they are 
either still far off shell, or on shell, but carry large components of the four-momenta, 
scaling like mq, or large energy release. The sum of the imaginary parts obtained 
in this way is the so called the quark-gluon cross section. This quantity serves as 
a reference quantity in formulating the duality relations. When one says that the 
hadron cross section is dual to the quark-gluon one, the latter must be calculated 
by virtue of the procedure described above. 

The contributions left aside in the above procedure are related, at least at a 
conceptual level, to the high-order tail of the power (condensate) series. They can be 
visualized with large-size instantons. A subtle point is that the large-size instantons 
contribute not only to the exponential terms, but also to the condensate (power) 
expansion. Our task is to single out the exponential contribution, since we have no 
intention to use the instanton model to imitate the low-order terms of the power 
expansion. Our instanton model is far too crude for that. In the next section 
we proceed to formulating the instanton model, making special emphasis on this 
particular element - isolating exponential contributions. 

Beyond the simplest one-variable problems, like the correlation function of two 
vector currents related to R e + e - or R T , very often one encounters a more complicated 
situation when the amplitudes have several separated kinematical cuts associated 
with physically different channels of the given amplitude, and one is interested only 
in one specific channel. This situation is typical for the inclusive heavy quark decays 



12, [L3f . The OPE-based predictions in this case require - additionally - a different 



type of duality: one needs to assume that a particular cut of interest in the hadronic 
amplitude is in one-to-one correspondence with the given quark-gluon cut. In other 
words, it is assumed that different channels (in terms of the hadronic processes and 
in terms of the quark-gluon processes) do not contaminate each other |Tj| . This was 
called "global duality" Q In the practical OPE the cuts of the perturbative coefficient 
functions carry clear identity, and the above assumption of "global duality" is easily 
implementable. The above assumption can actually be proven in the framework 
of the practical OPE at any finite order, as was shown in |13[ . However, most 
probably, this "global duality" fails at the level of the exponential terms. In the 
present paper we do not address the issue of "global duality" violations although 
instantons can model this phenomenon as well. Such effects are probably smaller 
than the deviations from local duality; in any case they deserve a dedicated analysis. 
It is worth noting that for one-variable problems (e.g. the total e + e~ annihilation 

1 Warning: the term global duality is often used in the literature in a different context. 
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cross section), duality for various integrals over the cross section over a finite energy 
range is still local duality, as will be discussed in Sect. 5. 



3 Abstracting General Aspects 

Before submerging into details of the instanton calculations, we outline the practical 
motivation for inclusion of the corresponding effects from the general perspective of 
the short distance expansion. It will also enable us to illustrate in a simple way 
the divergence of the power expansion. Consider a generic two-point correlation 
function I1(Q 2 ), say, the polarization operator for vector currents: 

n(Q) = J d 4 xe iqx Il(x) = - J d 4 xe iqx ( G(x, 0) G(0,x) ) (3) 

where G are the quark Green functions in an external gauge field and averaging over 
the field configurations is implied; we do not explicitly show the Lorentz indices. 
Equation (|3]), and all considerations in this section, refer to Euclidean space. 

The power expansion of n(Q) in 1/Q is the expansion of the correlation function 
in singularities near the origin. (This statement is not quite accurate in the 
Wilson's OPE; it is correct, however, in the practical OPE). Thus, one is interested 
in the small-x behavior of U(x) or, equivalently, of G(x). In the leading, deep Eu- 
clidean approximation, Green's functions are the free ones, ary/x 4 , plus perturbative 
corrections arranged in powers of a s (l/x 2 ): 

G pt (x) = JL ^ (i + ai a s (l/x 2 ) + ... ) ; x 7 = (4) 

(a particular invariant gauge is assumed here; the correlator is gauge-invariant any- 
way, and similar series can be written directly for the product of the two Green 
functions). All terms in the perturbative expansion (f|) have the same power of 
x, and differ only by powers of logx 2 . Logarithms emerge due to the singularity 
1/x 2 of the gluon interaction near x 2 = 0. Upon making the Fourier transform, the 
perturbative corrections in Eq. (Q) are converted into powers of logQ 2 . 

Power corrections, emerge from the expansion of Green's functions near 

x = 0: for example, in the Fock-Schwinger gauge 

where higher order terms contain higher powers of the gluon field, G a p, or its deriva- 
tives at x = (for a review see |I3|]). Since the additional terms in the expansion 



contain extra powers of x (generically accompanied by logx 2 ), it is clear that, re- 
turning to the momentum representation, one gets additional powers of Q in the 
denominator for extra powers of x in the expansion of U(x). (The positive powers 
of x in Il(x) are accompanied by logx 2 .) Thus, the 1/Q expansion obtained in the 
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practical OPE is in one-to-one correspondence, in the coordinate space, with the 
expansion of Green's function near the point x = 0. (This fact is absolutely explicit 
in ordinary quantum mechanics, where the dynamics are described by a potential. 
In QCD the power corrections to the inclusive heavy quark decay rates, for exam- 
ple, have a similar interpretation: the leading 1 / tuq correction, due to the Coulomb 
potential at the position of the heavy quark, is absent because of a cancellation 
between the initial binding energy and the similar charge interaction with the de- 
cay products in the final state. Physically, the reason is conservation of color flow. 
Moreover, the chromomagnetic term is determined by the magnetic interaction at 
the origin, and so on.) 

The question that naturally comes to one's mind is whether the above expansion 
in the x space is, in a sense, convergent. At best, it can have a finite radius of 
convergence which, for the given external field, is determined by the distance to 
the closest (apart from the origin) singularity in the complex x 2 plane. On the 
other hand, evaluating the Fourier transform (^) converting U(x) into H{Q), one 
performs the integration over all x. Therefore, even for arbitrarily large Q, one 
has to integrate U(x) in the region where the expansion of the Green's functions is 
divergent. Although any particular power term 1/Q n can be calculated and is finite, 
this leads to the factorial growth of the coefficients in the 1/Q expansion, and thus 
explains its asymptotic nature. 

Let us illustrate this purely mathematical fact in a simplified setting. Let us 
consider the "OPE expansion" of a modified Fourier transform (the one-dimensional 
integral runs from zero to infinity; such transforms are relevant in heavy quark 
theory, see ||): 

r°° 1 
f{Q) = / — e** = 

Jo x A + p z 

= y c ~ Qp + Y p [e ' QPWl{Qp) ~ eQPEi (-Qp)] . (6) 

where Ei is the exponential-integral function. The integrand has a singularity in the 
complex plane, at x — ±ip and is perfectly expandable at x — 0. Expanding the 
"propagator", l/[x 2 + p 2 ), in x 2 we get the "OPE series" 

/«) = /„ ^(-'^ = ? £^. < 7 > 

First of all note, that the "OPE series" has only odd powers of 1/Q. Comparing with 
the exact expression we see that the function f(Q) is not fully represented by 
its expansion (|7|), which is obviously asymptotic. The exponential term is missing. 
This exponential term comes from the finite-distance singularities of the integrand. 
Indeed, one can deform the contour of integration over x into the complex plane; the 
integral remains the same as long as the integration contour does not wind around 
the singularity at x = ±ip, whose contribution is 

5f(Q)=-e-^. (8) 
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This is precisely the uncertainty in defining the value of the asymptotic 1/Q series 
(0)- 

Since f(Q) is expanded only in odd powers of 1/Q, the symmetric combination 
g{Q) — f(Q) + f(~Q) has no power expansion at all. The function g(Q) does not 
vanish, however: 



/oo I p—Qp 
dz— — e«* = 7r . (9) 
-oo X + p p 



This expression is in agreement with the above estimate of the uncertainty of the 
power expansion per se and demonstrates that the exponential terms are present. 

The appearance of the terms exponential in Q in the example above bears some 
resemblance to the renormalon issue - the factorial growth of the coefficients in 
the perturbative expansions in QCD [T5 |. The Feynman graphs contain integration 



over all gluon momenta k 2 ; on the other hand, the expansion of a s (k 2 ) in terms 
of a s (Q 2 ) is convergent only for k 2 between some minimal and maximal scales, 
A.qcd ~ ^ 2 ~ Q 4 /A-qcd- Although each particular term in the expansion can be 
integrated from k 2 = to oo, yielding a finite number, the problem of divergence of 
the original expansion of a s (k 2 ) is resurrected as the factorial growth of the resulting 
coefficients. 

We conclude, therefore, that the divergence of the power expansion within the 
practical OPE, and the presence of the exponential terms is a rather general phe- 
nomenon, and is related, conceptually, to the singularities of Green's functions in 
the coordinate space at complex Euclidean values of x 2 located at finite distances 
from the origin. The question of the possible role of these finite distance singularities 
was first raised in Refs. JIB], |17| (see also [|T8| ). 



The analogy with renormalons in the a s perturbative expansions can be contin- 
ued. In the renormalon problem, the proper inclusion of the condensates, within the 
framework of the OPE, eliminates the infrared renormalons altogether and makes 
the infrared- related perturbative series well defined and, presumably, convergent. 
One may hope that a consistent explicit account for the finite-x 2 singularities would 
also make the infinite power series well defined. From the theoretical perspective, 
though, this problem has not been investigated so far. 

Addressing practical applications, there are two general reasons to expect that 
the inclusion of the exponential terms in the analysis can be important - even though 
the power series analysis accounts, at best, for only a few leading terms in the 
1/Q expansion. First, there exists some phenomeno logical evidence, to be discussed 
below, indicating that the impact of the exponential terms in the Minkowski domain 
may be more important numerically than that of the omitted condensate terms for 
intermediate values of the momentum transfers. This statement is illustrated by the 
t example, see Sect. 5. Second, historically, this is not the first case where we have 
encountered such a perverted hierarchy in QCD. It is quite typical that in the QCD 
sum rules, the contribution of the (omitted) higher order a s terms, which formally 
dominate over the condensates, is far less significant numerically than that of the 
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condensates. This point is crucial - while the exponential terms die out fast in the 
Euclidean domain, they decrease much more slowly in the physical cross sections 
and, thus, often dominate over the condensate effects, the more so that the latter 
often are concentrated at the end points, and are not seen at all outside the end point 
region. This observation was emphasized in Ref. [17]. Indeed, the terms ~ e~® p 
oscillate, rather than decrease, when analytically continued from the Euclidean to 
the Minkowski domain, Q — > iE. 

Leaving aside such subtle theoretical questions as the summation of the infinite 
condensate series, one may hope that including the principal singularities at the 
origin and at finite x 2 will lead to a description of the correlation functions at hand 
which is good numerically. Indeed, the leading singularity at x 2 = is given by the 
perturbative expansion, and subleading terms near x 2 = are given by the practical 
OPE. Adding the dominant singularity at x 2 ~ Aq^ d in the complex plane we 
capture enough information to describe the main properties of the function, and 
thus provide a suitable approximation to the exact result, which may work well 
enough for a wide range of x 2 , thus yielding the proper behavior of the correlation 
function, T1(Q), down to low enough Q 2 . One should clearly realize, however, that 
this procedure is justified only if we do not raise the subtle question formulated 
above and keep just a few of the first terms in both the perturbative and condensate 
expansions. Summing more and more condensate terms within this - rather eclectic 
- procedure may not only stop improving the accuracy, but even lead to double- 
counting of certain field configurations. Being fully aware of all deficiencies of this 
approach at present, we still accept it for estimating possible violations of local 
duality in a few cases of practical interest. Eventually, this approach may develop 
into a systematic, and self-consistent phenomenology of the exponential terms, much 
in the same way as the QCD sum rules represent a systematic phenomenology of 
the condensate terms. 

Technically, Green's functions are obtained by solving the equations of motion in 
the given background gauge field; in particular, quark Green's functions are obtained 
by solving the Dirac equation. Henceforth, the singularities outside the origin emerge 
only at such (complex) values of x 2 = where the gauge field is singular. Instantons 
provide an explicit example of such fields which are, of course, regular at real x 2 , 
but have a singularity in the complex x 2 plane 0. The singularities of the instanton 
fields are passed to the Green's functions, e.g. the singular terms in the spin-0 and 
spin-1/2 Green's function have the structure 

GJx, y) ~ -tit; , tit; , 

K ' {{x - z) 2 + p 2 f 2 ((y - z) 2 + p 2 ) 1 ' 2 ' 

Gf(x,y) 7 , , , £=-,-, (10) 

n UJ ((X-Z) 2 + P 2 f ((y_ Z )2 + p 2)* ' 2'2' 1 ' 

where z is the center of the instanton, and p is its radius. The actual nature of 
the singularity changes when one integrates over the position of the instanton, and 

2 The gauge potential itself may have singularities at real x 2 , but these are purely gauge artefacts. 
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over its orientation and size. The poles in U(x) may change into more complicated 
singularities (say, cuts). Interactions of different instantons will also affect the nature 
of the singularities. One does not expect, however, the finite-distance singularities 
to disappear. The precise position of the singularities, and their nature, depend on 
the details of the strong dynamics. 

The simplest finite x singularity in the physical correlator U(x — y) one can think 
of has the form 

U ( x ) = ( A 2Y ( H ) 
[x z + p l ) 

where v is some index. (The cases of v = 1 and 2 were discussed in Ref. fl7| .) Its 
momentum representation, 



n«« = / fe*n (l) = ^ (9iy- 2 !^m (12) 



clearly exhibits the exponential behavior related to the singularity at x 2 = —p 2 . The 
function on the right-hand side is exponentially small in the Euclidean domain but 
yields only an oscillating factor (damped by a modest power of 1/Q) upon analytic 
continuation to the physical domain, Q 2 = — s — iO, 



ImU(s) 



71-3 ( \fsp\ V 2 cos7rz/ J 2 -u(Vsp) + simru N 2 - u (y/ r sp) 



r» 2 p 



2^-4 



7T 3 ( y/spY 2 J u ~ 2 (Vsp) 
~ I» ^ 2 ) • ^) 

Note that at v = 1, the right-hand side of Eq. ([Tj^) implicitly contains S(s). Here, 
Bessel, McDonald, and Neumann functions are denoted by J, K, and N, respec- 
tively. Asymptotically, at large s, 



Im nW * -^:T4, (&) " "H^-(^^)f) . (14) 



5 5/4 y{v) \ 2 J p 

Certainly, the purely oscillating factor above is an extreme case. Any sensible smear- 
ing over p (with a smooth weight function) will restore the decrement of the exponent 
in the Minkowski domain (as discussed in Ref. ||). Generically, therefore, we obtain 
a decaying exponent, exp(— E a ), modulated by oscillations. More exactly, we get 
a sum of such terms. The index a depends on dynamics and, in principle, can be 
rather small numerically. If a is small, the damping regime takes over the oscillating 
regime at large values of E, after a few unsuppressed oscillations occur. At such 
values of E, the powers of 1/E in the pre-exponent can make the whole contribu- 
tion small. Therefore, starting our analysis with an extreme situation - a purely 
oscillating factor times some power of 1/E in the pre- factor - is quite meaningful. 
We will discuss all these details in a more specific setting of the instanton model. 
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It is easy to see that causality requires singularities of Il(a; 2 ) to lie either on the 
negative real axis of x 2 , or at larger arguments of x 2 , on the unphysical sheet; U(x 2 ) 
must be analytic at |argx 2 | < ir. In the instanton model the singularities are on 
the negative x 2 axis. Smearing the instanton sizes with a smooth function weakens 
the strength of the singularity near the purely imaginary y/x 2 and thus effectively 
moves it further into the complex plane. 

To summarize, we argued that the violations of local duality are conceptually re- 
lated to the divergence of the condensate expansion (practical OPE) in high orders. 
Technically, they may occur due to the singularities of Green's functions at complex 
Euclidean values of x at finite distances from the origin. Accounting for such singu- 
larities, in addition to the perturbative and the condensate expansion, is a natural 
first step beyond the framework of the practical OPE. In the next section we will 
proceed to a specific model for this phenomenon based on instantons. Since they 
are not necessarily the dominant vacuum component we try to limit our reliance on 
instantons to the absolute minimum. In particular, their topological properties are 
inessential for us, and even lead to certain superfluous complications. 



Here we will formulate our rules of the game. To get an idea of possible violations of 
duality we will consider a set of physically interesting processes (two-point functions 
of various currents built from light quarks, the transition operators relevant for the 
inclusive heavy quark decays and so on). Our primary goal is isolating the finite- 
distance singularities in x 2 , which will eventually be converted into the exponential 
terms in the momentum plane. To this end it will be assumed that the quark Green's 
functions in the amplitudes under consideration are Green's functions in the given 
one-instanton background. The one-instanton field is selected to represent coherent 
gluon field fluctuations for technical reasons - in this background Green's functions 
for the massless quarks are exactly known. 

The instanton field depends on the collective coordinates - its center, color space 
orientations, and its radius. Integration over all coordinates except the radius is 
trivial, and will be done automatically. Integration over the instanton radius p 
requires additional comments. 

First of all, in all expressions given below, integration over p is not indicated ex- 
plicitly unless stated otherwise. Any expression F(p) should be actually understood 
as follows 



where d(p) is a weight function and integration over instanton position, d^z/p 4 , is 
included in the definition of F(p). 

If we were building a dynamical model of the QCD vacuum based on instantons, 
we could have tried to calculate this weight function. As a matter of fact, for an 



4 Instanton Model 




(15) 
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isolated instanton the instanton density, d(p), was found in the pioneering work [19 
for pure gluodynamics, 

d(p) Q = const (pAqcd) 6 , (16) 

where b is the first coefficient in the Gell-Mann-Low function (b = 11/3 iV c for the 
SU(3) gauge group). Of course, the approximation of the instanton gas |2(| is totally 
inadequate for many reasons - one of them is the fact that inclusion of the massless 
quarks completely suppresses the isolated instantons [[HJ. This particular drawback 
can be eliminated if one takes into account the quark condensate, (qq) ^ 0. Then 
the instanton density takes the form plj 

d(p) = const ((qq)p 3 rd (p) (17) 

where n/ is the number of the massless quarks, and now b = 11/3N C — 2/3nf. Note 
the extremely steep p dependence of the instanton density at small p. The impact 
of the quark condensates is not the end of the story, however, since for physically 
interesting values of p, the vacuum field fluctuations form a rather dense medium 
where each instanton feels the presence of all other fluctuations. In principle, one 
could try to build a model of the QCD vacuum in this way, for instance, that is what 
is done in the so called instanton liquid model (see f22| , |23| and references therein). 
The main idea is that the instanton density is sharply peaked at p ~ 1.6 (GeV) -1 , 
where the classical action is still large, i.e. we can still consider individual instantons. 
On the other hand, the interaction between instantons is also large, but still not large 
enough so that the instantons melt. The extremely steep growth of the instanton 
density at small p is cut off abruptly at larger p, due to interactions in the instanton 
liquid. The proposed model density which captures these features is just a plateau 
at p = p c with the width 5 <^ p c . 

We would like to avoid addressing dynamical issues of the QCD vacuum in the 
present paper. Our task is to rely on general features, rather than on specific details, 
and the instanton field, for us, is merely representative of a strong coherent field 
fluctuation. For this limited purpose, we can ignore the problems of the calculation 
of the instanton density, and just postulate the weight function d(p) in the simplest 
form possible. The most extreme assumption is to approximate d(p) by a delta 
function, 

d(p) = d p 5(p - p ), (18) 

where do and po are appropriately chosen constants. In a very crude approxima- 
tion this weight function is suitable, in principle, although it has an obvious draw- 
back. If p is fixed, as in Eq. (|18|), the instanton exponential in the Euclidean 
domain becomes cosine in the Minkowski domain, with no decrement. For instance, 
(Qp)~ 1 Ki(Qp) — > E~ z l 2 cos (Ep — phase), where the arrow denotes continuing to 
the Minkowski domain, taking the imaginary part, and keeping the leading term in 
the expansion for large Ep. If one wants to be more realistic, one should introduce 
a finite width. A reasonable choice might be 

«;(p) = (p)-Wexp{-- -/?/>} (19) 

P 
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where TV is a normalization constant, 



p is the center of the distribution, and A is its width. Convoluting (Qp)~ 1 Ki(Qp) 
with this weight function, one smears the cosine, which results in the exponential 
fall off in the Minkowski domain, 



(Qp^K^Qp) -> AT2E- 1 MV^l^P 2 + £ 2 - flWyidiVtetyP 2 + £ 2 + 3} 1 ' 2 } 

(21) 

where the meaning of the arrow is the same as above. 
If 

E » i^, (22) 

the imaginary part reduces to 



2NE-\h{pEp7f)K l {^2Ep~, P f) , (23) 
and falls off exponentially. If the weight function is narrow, (A -C po), this exponen- 



tial suppression starts at high energies, see Eq. (22). In the limit when A — > 0, with 



E fixed, the exponential suppression disappears from Eq. (2TJ), and we return to 
the original oscillating imaginary part. Note also that the exponent at E ^> is 
different from the one in the Euclidean domain {\fE versus Q). In Sect. 5.2 we will 
introduce the corresponding index, a, characterizing the degree of the exponential 
fall off in the Minkowski domain at asymptotically large energies. 

Concluding this section, we pause here to make two remarks of general character. 
The fact that smearing the scale with smooth functions of the type (|T9|) produces 
exponential fall off is not specific to the instanton-induced spectral density. Even 
much rougher spectral densities (with appropriate properties), being smeared with 
the weight function ((03), become exponential. An instructive example is provided by 
a model spectral density suggested in Ref. || . Consider the following "polarization 
operator" 

™(^) 

where 3 is the special beta function related to Euler's ip function, 



m = \ 



^TTk- (24) 



This fake polarization operator mimics, in very gross features, say, the difference 
between the vector- vector and axial-axial two-point functions |. At positive Q 2 , it is 
expandable in an asymptotic series in 1/Q 2 , plus exponential terms. At negative Q 2 , 

3 In Ref. [|| the model spectral density (|24|) was suggested in the context of the heavy-light 
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(positive s), it develops an imaginary part. The imaginary part obviously consists 
of two infinite combs of equidistant delta functions - half of them enter with the 
coefficient 1, the other half with the coefficient —1. Literally speaking, there is no 
local (point-by-point) duality at any energies. 

Let us smear the combs of the delta functions with the weight function (|19|). 
Now, the imaginary part at negative x is smooth, exponentially suppressed, and 
oscillating, 



fOC <jr ^ j / 

Im / dp w(p) (3{px) = -Y(-l) k w(k/\x\) oc Ime-V-/ 2 MV«N+^. (25) 
Jo x ^ 

Indeed, one can represent the sign alternating sum, (— l) k w(k/\x\), as the integral 
of the function i/(2sm(iTz) w(z/\x\)), with complex variable z, over the contour 
embedding the positive real axis [l,+oo). At large \x\, its value is determined by 

large z; the integrand has two complex conjugated saddle points, z = \f^^ e ±l7r/ ' 4 , 
whose steepest descents lead to z = and z = ±ioo. Evaluating the saddle point 
integrals, one arrives at the above asymptotics. 

Returning to the instanton model, we note that the weight function, (0), is con- 
venient, because the contribution of the small-size instantons (which affect the OPE 
coefficients and are not discussed in the present paper) are naturally suppressed. 
The absence of these small-size instantons allows for a sensible expansion param- 
eter, l/(mgp), which can be used in calculations with heavy quarks. It is worth 
emphasizing again that at very large momentum transfers (energies, heavy quark 
masses, etc.), the small-size instantons will always dominate over the exponential 
terms. Thus, our model is applicable, if at all, only to intermediate scales. 

In QCD, the instanton field configuration does not constitute any closed ap- 
proximation. Therefore, one may question practically every aspect of the model we 
suggest. Developing phenomenology of the exponential terms will help us under- 
stand whether this approach has grounds. From the purely theoretical standpoint it 
might be instructive to consider a formulation of the problem where the instantons 
can be studied in a clean environment, rather than in the complex world of QCD. 
Such an analysis was already outlined in the literature [25]. Let us assume that in- 



stead of QCD, we study the Higgs phase, i.e. we introduce scalar colored fields which 
develop a vacuum expectation value, and break color symmetry spontaneously. The 
gluon fields acquire masses. If their masses are much larger than Aqcd, we are in the 
weak coupling regime, and the semiclassical approximation becomes fully justified. 

two-point functions, and, accordingly, x was related to E, not Q 2 . In the heavy-light systems, 
the model does not reproduce fine features either; in particular, the equidistant spectrum it yields 
is not realistic. The separation between the highly excited states should fall off as l/E. Such a 
behavior immediately follows from the semiclassical quantization condition 



(E — A r)dr oc n . 

Previously this pattern was noted in the two-dimensional 't Hooft model 124 
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The instanton contribution to various amplitudes is well-defined now, and subtle 
questions, which could not be reliably answered in QCD, can be addressed. 

The pattern of the instanton contribution as a function of energy in this case was 



studied in Ref. |po] . It is quite remarkable that the pattern obtained bears a close 
resemblance to what we have in QCD, in particular, oscillations in the Minkowski 
domain. 



5 Deviations from duality in R e + e ~ or hadronic r 
decays 

5.1 Instanton estimates 

The peculiar details of local duality violations are more transparent in the simple 
cases of e + e~ annihilation cross section and inclusive hadronic r decays. Several 
rather sophisticated analyses of the instanton effects in these problems were carried 
out recently |26], |27|, |28|, E9| . A more general consideration, rather close in ideology 



to our approach, was given in [jTjJ (in a sense the spirit of the suggestion of Ref. [17 
is more extreme). We further comment on these works in Sect. 11. To see typical 
features of the instanton-like effects we consider, for simplicity, the correlator of the 
flavor-nonsinglet vector currents relevant to R e + e -; a similar correlation function 
appears in the hadronic r decays, alongside with its axial-vector counterpart. For 
simplicity, we will mainly ignore the latter contribution and discuss the vector part 
as a concrete example. 

Let us define (Q 2 = —q 2 — iO) 

IV(g 2 ) = / rf 4 xe^(zT{j;(x)J,(0)}) = ^ ( q \ u - q ^ u )U(q 2 ) , 

U(Q 2 ) = log% + ... , (26) 
I* 1 

R(s) = --Imn(Q 2 = -s -iO) = 1 + .... 

7T 

For the purpose of our discussion the average (...) is not yet understood as averaging 
over the physical vacuum; we, rather, calculate the correlation function in a partic- 
ular external field, and average over certain parameters of this field (the invariant 
tensor decomposition is appropriate in the latter case). The second equation shows 
Tl(Q 2 ) in the absence of any field. In a given field, IT^ (x, y) is merely a trace of the 
product of the two Green functions which are explicitly known for massless quarks 
in the field of one instanton. Upon averaging over the positions and orientations of 
the instanton of the fixed size p, one arrives at the known expression (integration 
over p is shown explicitly) 



n(Q) = n (Q) + u\Q) 
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1 ^ 2 

fl 2 



+ 16tT 



dp 



d{p) 



dtKo 



2Qp 



(27) 



.3(Qp) 4 (QpfJo —\^T—t2 / 

where K 2 is a McDonald function, and J(d(p)/p 5 )dp is to be identified with the 
number of instantons per unit volume. The superscript I marks the instanton con- 
tribution. The first term in the square brackets is "a condensate" , the second one, 
on the contrary, does not produce any 1 / Q n expansion. Considering Eq. Q27D in the 
Minkowski domain one has (E = y/s) 



R{E) = Rq(E) + R\E) 



■57T 



dp 
P 



d{p) 



-AE 2 ) + 



dt Jo 



2Ep 



t 2 



(28) 



2p 2 ""~ ' ' (Ep) 2 J VVT 

where J2 is a Bessel function. Violation of local duality at finite E is given by the 
last term. Expanding the Bessel function at large Ep, and performing the saddle 
point evaluation of the inner integral, we see that it oscillates, but decreases in 
magnitude only as 1/E 3 : 



R\E)\ 



fixed p 



-An 4 



1 



(E P y 



cos (2Ep) 



(29) 



(we remind the reader that the true power corrections from the OPE appear in the 
imaginary part at large E only at the level ce^/E 4 provided that the quarks are 
massless, as we assume here). After averaging this result over p with a smooth 
enough weight, the resulting R 1 decreases exponentially at E — > 00; the decrement 
is determined by the analytic properties of d(p). The behavior of the "exponential" 
contribution, given by the last term in Eq. (p8|) at small E, is relatively smooth, 



R\E) 4tt 2 log (Ep) 



(30) 



To visualize the above expressions, we plot in Fig.|l| the value of R(E) stemming 
from Eq. (|28| ) in the real scale E using the instanton density (0), with some rather 
ad hoc overall normalization f] do, and po = 1.15 GeV -1 . 

represents experimental values extracted from the CLEO data |3~T 



Figure 



Although our theoretical curve does not literally coincide with the actual data, 
it is definite that the general feature of the experimental curve - the presence of 
oscillations, and a moderate falloff of their magnitude with energy - is captured 
correctly. The fact that our model is not accurate enough to ensure the point-by- 
point coincidence was to be anticipated. Obvious deficiencies of the model will be 
discussed in Sect. 10. Some additional remarks concerning duality violations in the 
hadronic r decays are given in Sect. 9.1, see Eq. ( |105| ). 

4 In Sect. 9.1 we discuss our choice for do. Our motivation is based on phenomenological analysis 
of the duality violations in the semileptonic D decays, which may be as large as ~ 50% . 
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Figure 1: R(E), taking into account the instanton contribution. The perturbative 
result is normalized to unity. 




Figure 2: Experimental value of R(E). 
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5.2 Three zones 



A single glance at experimental data (a part of the data is presented on Fig. 2) 
reveals a striking regularity of the inclusive cross section. We believe that this 
regularity is a general phenomenon, and its discussion is very pertinent to the issue 
of the duality violations. The same pattern of behavior is expected, say, in the 
spectra of the radiative decays B — > X s + 7, and so on. 

One can single out three distinct zones in the physical inclusive hadronic cross 
sections, governed by different dynamical regimes. If we proceed from the low in- 
variant masses of the inclusive hadronic state to high masses, the first zone we see is 
a "narrow resonance" zone. It includes one, or at most two, conspicuous resonances. 
It stretches up to a first boundary - call it s . Crossing this first boundary, we find 
ourselves in the second zone - the oscillation zone. The cross section here is already 
smooth, and the point-by-point violations of the quark (gluon) hadron duality are 
not violent. Still, these violations are quite noticeable (they may constitute a few 
dozen percent), and have a very clear pattern - several clearly visible oscillations, 
with relatively mild suppression, 



The upper boundary of this zone will be referred to as si. Finally, above this second 
boundary, there lies a third domain - the asymptotic zone, where 



R = i? O PE + exp[-(2£p) CT ]sin((2£p , r + 0) or {l/E)\ a < 1 , 7 >1. 



Here Rope is a smooth (practical) OPE prediction, k is an integer, a and 7 are 
indices. Our model is intended for applications in the second (oscillation) zone. 

It is worth noting that the precise values of the boundaries sq and s\ are very 
sensitive to dynamical details. For instance, in the imaginary world with infinite 
number of colors, s , is believed to go to infinity, and the regime of the second zone 
never occurs. 

5.3 Smearing and local duality. 

In this section, we discuss another general, and crucial feature of the "exponential" 
terms. What happens if, instead of considering the imaginary parts point-by-point, 
we choose to analyze some integrals over a finite energy interval, with some weight? 
Intuitively, it is clear that violations of the quark (gluon) - hadron duality are 
expected to become smaller if the weight function is smooth enough and the energy 
interval over which we integrate is large. The case when one integrates with a 
polynomial weight (polynomial in s = E 2 ) is of a particular practical interest. Let 
us consider the finite-energy moments 



R = R OPE + (const/£ fc ) sin(2£p + 0) . 




(31) 
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The deviations from duality are smallest at the upper edge of the integration 
domain, and largest at the lower edge. Intuitively, it is clear that the deviations 
from duality in the integral (|31|) are determined by deviations at the upper edge 
of the integration domain. This result, however, is obtained only if one makes full 
use of the analytic properties of the exponential contributions at hand. If one tries 
to directly integrate the asymptotic instanton formulae over t, in a straightforward 
manner, one gets a huge contribution determined by the lower end. This is the 
essence of the so called "a part larger than the whole" paradox, observed in the 
instanton calculations, say, in Refs. J7|, ||, where the instanton contribution to the 
decay spectrum turned out to be parametrically larger than that to the total decay 
rate. 

Let us elucidate the point in more detail. The moments, JA n (s), get contribu- 
tions both from the usual OPE terms, which are located at small s ~ ^-qcd ( m 
our case it is 5'(s) from the term 1/Q 4 in Eq. (p7|) which survives only for n = 1), 
and from the exponential part going beyond the practical OPE. We are interested 
here only in the latter piece and, therefore, subtract the condensate part. Using the 
large-s expansion for R, one is literally in trouble: the integral over the imaginary 
part ( p9|) seemingly diverges at small s where this expression is not applicable, and 
must be cut off at s < 1/p 2 - At first sight, it then seems the result completely 
depends on the lower limit sq, and on the precise way of implementing the cut off 
at So- The fact that we integrate over a large interval stretching up to s seems to 
be of no help in suppressing the duality violations. 

It is easy to see, however, that the large result above is obtained only because 
we have used a wrong expression at small s. The asymptotic instanton formula is 
definitely invalid at small s. If p is fixed, we could use, of course, the exact instanton 
expression at small s which is (almost) not singular (see Eqs. ( |3~0l) and (|28|) ). We 
would not trust the instanton result at small s anyway. Therefore, the prediction for 
the moments, Ai n (s), should be obtained without relying on the expicit expressions 
at small s. To this end one invokes dispersion relations. 

Whatever the origin of the exponential contribution under consideration is, it 
must obey the dispersion relations. Take H(Q) — Hope{Q), where the "practical 
OPE" piece, IIope(Q), in our example is given explicitly by the single term 

(32) 

(we omit the subscript I, since, in what follows we consider only instanton induced 
contributions). Since U(Q) — IIope(<3) exponentially decreases at large Euclidean 
Q 2 , one has an infinite number of constraints 

lim Q 2n (n(Q) - IWQ)) = (-1)" / (R(s) - R OP e(s)) s^ 1 ds = . 

(33) 

In other words, all moments considered in the full s range from to oo, Ai n (oo), are 
given completely by their OPE values, and the extra contribution from the non-dual 
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piece is absent. (To define A4 n = A4 n (oo) in the particular example one may need 
to regularize integrals in (ftTD, (|33| ) by, say introducing a damping exponent e~ €v/i 
with an infinitesimal e.) 

Using this property, one immediately concludes that the violation of duality in 
the moments, A4 n (s), is determined, parametrically, by the upper limit of integration 

^ ' POO 

M n (s) = M° PE (s) - (n + 1) / R(t)t n dt . (34) 



This is, clearly, the most general property of the "exponential" terms which does 
not depend on any details of a particular ansatz. 

Formally, the relations of the type (|33|) and (|34]) for the imaginary part (obtained 
by the analytic continuation of the Euclidean exponential terms) can be written as 
follows: _ 

r r 81 

R(s)-R pe(s) = / dt$(t) 5(t-s)-e- t a-s5(s) , (35) 
Jo 1 J 

where $(£) vanishes at t < and coincides with the asymptotic instanton expression 
at positive t, 

$( s ) = -4vr 2 -J_ cos (2 Pv ^) + O (V 2 ) . (36) 

In other words, to do the smearing integrals properly one must substitute R(s) by 
R(s) plus the whole tower of terms presented on the right-hand side of Eq. (j35|). 

The representation (^) is convenient since it explicitly ensures the property (0), 
which is the fact that the corresponding contribution to the correlator dies out faster 
than any power in the deep Euclidean domain. It shows that any particular spectral 
density generated at large s as a violation of local duality, must be necessarily accom- 
panied by the corresponding OPE-looking terms located at small s; disconnecting 
these seemingly different contributions is not consistent with analyticity. 

Let us parenthetically note that similar relations, with delta functions at the 
end point, must be used in the instanton calculations of the semileptonic spectra 
in the heavy quark decays. The occurrence of the end-point delta functions in the 
instanton expressions is reminiscent of what happens with the regular (OPE) power 
corrections to the semileptonic widths [0. The interaction with the final quark 



magnetic moment does contribute to the inclusive lepton spectrum with a definite 
sign in its regular part. And, yet, it is known to be absent in the total width. The 
cancellation occurs due to the terms located at the end point of the spectrum which 
- in the naive approach - are not seen in the 1/rriQ expansion. 

Loosely speaking, the part referring to low t in the integral (|3l|) is eaten up by 
the "condensates". 

Eq. (^) demonstrates that the deviations from duality in the finite-energy mo- 
ments, A4 n (s), are generically given by the accuracy of local duality at the maximal 
energy scale covered, s. More exactly, the error is approximately given by the inte- 
gral over the last half-period of oscillations. For the sign-alternating combination of 
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the moments, similar to the one determining the hadronic width of the r lepton, 



r had (r) ~ 2M Q (m 2 T )/m 2 T - 2M 2 (m 2 T )/m e T + M 3 (m 2 T )/m 8 T , 

it is likely to be larger and can be governed by a lower scale. Indeed, the resulting 
weight function 

2 > /- - s 




w T (s) = 2 #(m 2 T - s) • 1 + 2— 1 - — (37) 



is saturated mainly at s < m 2 /3. Therefore, s ~ m^/3 can be viewed as the actual 
mass scale governing the duality violation in this problem. Numerically it is close 
to the first pronounced resonance in the axial channel. 

An obvious reservation is in order here. In real QCD, where the series of the 
power corrections in the OPE is infinite, and presumably factorially divergent, un- 
tangling the exponential terms from the high-order tail of the series remains obscure. 
There is no answer to the question "what is the summed infinite OPE series?" , even 
in the Euclidean domain. Our approach to this issue is purely operational, and is 
clearly formulated in simple problems: pick up the contribution of the finite x gular- 
ity in the saddle point approximation. It is motivated by the general consideration 
of Sect. 3. 



6 Soft instantons in the 1/rriQ expansion. 

In this section, we briefly outline the generalities of the instanton induced expo- 
nential corrections of heavy quark decays. The goal of this section is a "back of 
the envelope" calculation presenting the functional dependence on the heavy quark 
mass. More detailed calculations, which will provide us with all coefficients in the 
pre-exponential factors, are deferred until Sects. 7 and 8. 

The main feature of the problems at hand is the presence of a large parameter, 
rriQp, which allows us to obtain sensible analytic expressions. As we have already 
discussed in a general context, there are three types of contributions associated with 
instantons: (i) Small size instantons affect the coefficient functions of the OPE. We 
are not interested in these terms. They will not appear in our calculations, since 
instantons of small size are, by definition, excluded from our model density function 
d(p), and we always assume that rriQp ^> 1. Technically, as was already mentioned, 
small size instantons cannot be taken into account using the standard methods 
of HQET. (ii) The terms proportional to powers of 1/(toqp). They represent the 
instanton contributions to the matrix elements of various finite-dimension operators 
that are present in the OPE. (In the present context these terms are actually pure 
contamination, and so we will discuss only how to get rid of them), (iii) Finally, 
there are exponential terms, of the form exp (— 2rriQp). These terms are our focus. 
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Figure 3: Forward scattering amplitude. Bold lines represent propagation of a 
particle in the instanton field. 

6.1 Decays into light quarks 

Consider the generic form of an inclusive forward scattering amplitude which corre- 
sponds to the decay of a heavy quark into a massless quark and a number of color 
singlet particles, say, Iv (Fig.|), 



Here x and y are the position of the heavy quarks, and z is the instanton center (the 
integration over x + y yields the (^-function in the transition amplitude expressing 
the conservation of the total 4-momentum, which we do not write explicitly). Q(x) 
is the field of a heavy quark with mass itiq in the instanton background, S(x, y) is 
the Green function of the massless quark in the instanton background, and G s (x — y) 
represents the product of all color singlet particle (non hadrons) Green's functions 
produced in the decay. Note that all Lorentz indices are suppressed, as well as the 
integration over the instanton parameters, other than its center. For simplicity, we 
do not explicitly indicate the dependence of the fields and the quark Green function 
on the instanton collective coordinates, except for position. In the following, we will 
use the singular gauge for the instanton field. In principle, speaking of instantons 
assumes that expressions are written in Euclidean space, but so far the exact nature 
of the external field is inessential. 

We will always assume the heavy hadron is at rest, and thus we can single out 
the large "mechanical" part of the x-dependence in Q(x) : 



Calculating the width we will need to calculate the expectation value of the transition 
operator between the heavy hadron state: 



T = i J Q(x)S(x,y)Q(y)G s (x - y)d 4 (x - y)d A z 



(38) 



Q{x) = e - im Q l 



Q(x) 



(39) 



T 



1 



(H Q \T\H Q ) 



(40) 



2M, 



where now 



f = i jQ(x)S(x,y)Q(y)G s {x-y)e im ^ x °- yo) d 4 (x-y)d 4 z . 



(41) 
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Since the product of the quark Green functions and nonrelativistic Q fields does not 
have an explicit strong dependence on itlq, we clearly deal with a hard (momentum 
~ tuq) Fourier transform of a certain hadronic correlator which is soft in what 
concerns nonperturbative effects. Note that rriQ can now be considered an external 
parameter in the problem, for example, as an arbitrary, and even complex, number. 
We are not yet formally ready, however, to consider the Euclidean theory since we 
still have initial and final states. We shall address this issue a bit later, and now 
proceed as if we deal with free heavy quarks, which are transferred to the Euclidean 
domain without problems. 

Let us examine the propagator of the massless quark in the instanton back- 
ground, which is calculated exactly for the case of spin 0, 1/2, 1 particles [|33|. This 
(Euclidean) Green function S(x, y) has a generic form 

S ( X > V "> = [(a; - y)2]n [( x _ z )2 + p 2]fc l/2 [( y _ z )2 + p 2]fc 2 /2 X ^ (42) 

where S has no singularities at complex x a or y a (a polynomial). Using the Feynman 
parametrization, we rewrite it as 

Q( r y ) = r ( fc ) 1 5 x 

* [X ' y) ^ IW2)IW2) [(x-y)r 

I" dU kl/2 -\l-i) k *' 2 - 1 TTTT^v WtTT ^ > ( 43 ) 

Jo [£(1 - C) (x - y) 2 + (p 2 + z 2 )f 

k= h + k 2 ^ z = z-£x-(l-€)y . 

For the propagator of a spin or 1/2 particle, the value of k is eventually 1 and 2, 
respectively, and n = 1, 2. The large-momentum behavior of the Fourier transform 
of the correlator, Eq. (|4T|) , depends on the analytic properties of the integrated 
function. Let us first consider the analytic properties of S(x, y) in the complex 
(xo — yo) plane (Fig.^) - it has two different singularities. One singularity is on the 
real axis, and corresponds to two quarks being at the same point. This is the same 
singularity occuring in the Green function of free quarks, but upon integration, the 
residue is softly modified by the instanton field. Picking up this pole and calculating 
the amplitude, we will get instanton contributions to the usual power (l/mgp) n 
terms in the OPE, which we are not interested in. Indeed, making a Taylor expansion 
in (x — y)/p around this pole, we obtain a series of corrections ((x — y) / p) k / (x — y) 2n , 
which, integrated with the exponent, result in the above terms. 

Another singularity lies on the imaginary (xq — yo) axis. It comes from the finite 
quark separation 

(x-y) 2 =m-0}-\p 2 + ^) ■ (44) 

In contrast to the perturbative or OPE pieces, this separation does not scale like 
1/ttiq, but stays finite in the heavy quark limit. Upon integration over d i zd i {x — 



25 



Im x 

C) iP 



Re x 



Figure 4: Finite distance singularity of the Green function of a massless particle in 
the instanton background. 

y) <i£ this singularity, together with the factor e im Qi x o-yo) f r0 m the heavy fields, 
produces the e _const m Qf terms in the (Euclidean) amplitude that we are looking 
for. We then only need to determine the constant that enters the exponent, and the 
pre-exponential factor. 

Now with this general strategy in mind, let us outline the machinery in more 
detail. We want to abstract from the complicated questions of the interrelation 
of the instanton configurations to the particular heavy hadron structure, i.e. to 
consider the simplest possible state similar to a quasifree heavy quark instead of a 
real B or D meson or heavy baryon. On the other hand, the heavy quark, a priori, 
cannot be taken as free since the field Q(x) must obey the QCD equation of motion, 
in particular, in the instanton field. It is clear that such a program can be carried 
out consistently if the instanton size is small enough compared to the typical size of 
the hadron ~ Aqcd, but still is much larger than ttiq 1 . Having this choice in mind, 
we neglect, in what follows, the fact that the heavy quark is actually bound in the 
hadron although, eventually, the values of p will not be parametrically smaller than 
the hadronic scale. 

Thus we merely solve the equation of motion for the heavy quark in the instanton 
field, as one would do for an isolated particle. The role of the initial hadronic state, 
Hq, in Eq. (f|0j) is played by the single heavy quark spread in space and evolving 
in time according to the solution of the heavy quark Dirac equation analytically 
continued from Euclidean to Minkowski space. In our actual calculations we, of 
course, go in reverse: both the heavy quark field and the transition operator are 
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calculated in the Euclidean domain; the subsequent continuation to the Minkowski 
space is performed in the final expression for the forward amplitude T. Technically, 
we are able to solve the equation of motion for the heavy quark field since the 
parameter rriQp 1. 

The heavy field Q(xq, x) can be written in the leading order as 

Q(x ,x) ^Te'CW Q(0,x)+O(l/(m Q p)) = 

U(x)Q(0,x) + O(l/(m Q p)) . (45) 

The expression is written in Euclidean space, although we use Minkowski notations. 
Using the explicit solution for the SU(2) instanton in the singular gauge, Eq. (|62]), 
one gets the matrices U in the following cumbersome form: 

U(x) = exp lirft 




where ft = (x — z)/^(x — z) 2 , z is the coordinate of the instanton's center and 
t/2 are the color 577(2) generators. The integration is simplified since along the 
integration path x — z = const in Eq. fl4*5|) , Aq is proportional to one and the same 
color matrix f{x — z ) and, therefore, the path-ordered exponent reduces to the usual 
exponent of the integral of Aq. 

It is important that the expression for Q(x) is only valid in the leading order 
in the expansion parameter l/mqp. If one considers the contribution of small size 
instantons (as in Ref. 0), then no legitimate expansion parameter is available. The 
expansion in the heavy quark mass can only be obtained if instantons of size p < p c 
(where p c is some parameter ^> l/mq) are absent. Otherwise, the corresponding 
equations of motion need to be solved exactly. 

Even though we managed to solve the equations of motion for the heavy quarks 
in the leading approximation, the solution governed by the color matrices U is not 
analytic. The apparent singularities at x = z or y = z are, in fact, spurious and 
merely an artifact of using the singular gauge for the instanton field; since the 
amplitude we calculate is manifestly gauge invariant (it is nothing but the light 
quark Green functions times the path exponent []), this singularity is absent in the 
full expression, being canceled by similar terms in the light quark Green functions. 

5 Since the path exponent over the straight line is unity in the Fock-Schwinger gauge, the 
products we calculate are nothing else than the quark Green functions in the instanton field in the 
Fock- Schwinger gauge. However, the fixed point of the gauge does not lie at the center of the 
instanton, but, rather, at the external current point, x or y. For a review of the Schwinger gauge 
see Ref. Q. 



arctan ( 7—^—7 I - arctan 
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However, in general, the propagation matrix U introduces additional exponential 
corrections, since 

arctan ( ■== ) 
\ \Jx l + p 2 j 

has a (cut) singularity at t = i^x 2 + p 2 , which is a point where the singularity in 
the potential is not of the gauge type. The fortunate simplification which arises, 
in the leading approximation in l/(mQp), is that the two factors U(y) and U^ 1 {x) 
are unity at the saddle point. This happens due to the fact that the saddle point 
corresponds to the configuration where the instanton is situated right on the line (in 
three-dimensional coordinate space) connecting x and y and thus x — z = y — z = 
Eq. ([48]) Finally, picking up the pole in the complex (xo — yo) plane in the light quark 
propagator in Eq. ( ]45| ) , 

[((*-y) 2 + K(i-orv + 5 2 )]- fc , (47) 

we get an exponential factor 

exp (-m Q j{p> + F)/(£(l - 0) + (x- y) 2 ) 

in the transition amplitude. The expression in the exponent has a sharp minimum 
at 

5 = 0, £ = 1/2, (x-y) 2 = . (48) 
Evaluating it at this point we get the exponential factor 

e -2m Q p 

The power of ijiq in the pre exponent can be determined without actual calcula- 
tions as well. The residue of the fc-th order pole of the propagator yields the factor 
m O~ 1 upon integration over x — y . The Gaussian integrals over (x — y), xq — yo, z, 
and £ around their saddle points give ttiq 4 ; on the contrary, all "free" propagators 
(those of the color singlet final particles and the bare propagators of the quarks 
produced) enter at a fixed separation ~ — p 2 , and are mg-independent. We thus get 

T oc const • m k Q - 5 e ~ 2m ^ p . (49) 

For example, the large size instanton corrections to the transition amplitude for the 
semileptonic decays of heavy quarks has the form 

e -2m Q p 

T A oc 3— . (50) 

m Q 

The pre exponent can be easily calculated in the same way and will be given for a 
few cases of interest in the subsequent sections. 
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Figure 5: Exact (solid line) and asymptotic (dashed line) behavior of R(E). 



The same counting rules apply for the case when more quarks are present in 
the final state. In the case of the vector correlator of the light quarks we have the 
product of two light quark propagators instead of one, and eventually k = 4. Since 
IWQ) ~ Q 2 n(Q 2 ), we get for IT(Q 2 ) defined in Eq. fl27|) 

nw2> = 4 * 3 fiS ■ <5i) 

in accordance with the explicit calculation of Eq . (p7|) . In fact, the asymptotic 
expression works accurately enough even in the Minkowski domain already at y/sp ~ 
3; the corresponding approximate expression for -Imll is plotted as a dashed curve 
in Fig.|5|. It is clear that the main effect of the subleading in l/(Qp) terms is a phase 
shift in the oscillations. 

Using the same counting techniques, one concludes that the non-OPE soft in- 
stanton contribution to the forward amplitude describing nonleptonic decays of the 
heavy quark (we have three light quarks in the final state i.e. k = 6) scales like 

T nl ex rnqe"^' . (52) 

In the next section, we give a more detailed calculation, and discuss total widths 
and differential distributions. 

A qualifying comment is in order here. So far, in calculating the exponential 
terms, we assumed the heavy initial quark to be static (at rest). In this approxi- 
mation the presence of the initial heavy quark does not affect our result at all; it 
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is the final quarks that fully determine the exponential terms. We know for sure, 
however, that the initial heavy quark experiences a "Fermi" motion inside the heavy 
hadron. In the practical OPE, the first correction due to this Fermi motion comes 
from the operator Qtx 2 Q. In other words, the accuracy of the static initial quark 
approximation is 1/rriQ. An important question is how the exponential terms of the 
type we focus on show up in those subleading effects proportional to (n 2 ) (and other 
similar subleading effects). 

The expansion of the heavy quark propagator in I/tuq leads, generally speaking, 
to terms of the type (l/mg)l/(a; 2 + p 2 ), i.e. we get an extra singularity in the 
amplitude due to the initial heavy quark. This extra factor enhances the overall 
singularity of the amplitude in the complex plane and, thus, leads to a higher power 
of rriQ in the preexponent. The additional denominator, (x 2 + p 2 ), generates only the 
first power of rriQ, however, so the overall effect is suppressed by the small parameter 
l/(m Q p). 

In the Minkowski domain, inside the oscillation zone, the exponential factor is 
not a suppression at all, so we must count only the pre-exponential factors. We see 
that deviations from duality are parametrically relatively stronger in the (-7T 2 ) piece. 
Still the original suppression of the (n 2 ) piece by 1/toq is not completely lifted. We 
lose one power of 1/rriQ, but retain the other power. Thus, our approximation of 
the static initial quark in the analysis of the exponential deviations from duality is 
justified. Nevertheless, it is interesting to note that the initial-state I/tuq effects 
are less suppressed in the exponential terms. This seems to be a general feature. 

6.2 Decays into massive quarks 

The case of the massive final state quark (e.g. b — > civ) does not differ conceptually 
if treated in the l/(Qp) expansion, although a technical complication arises due to 
the unknown explicit expression for the massive propagator in the instanton field. 
Still, this is not a stumbling block in the analysis: the relevant singularity of the 
c quark Green function can again only be at (x — z) 2 = —p 2 or (y — z) 2 = —p 2 , 
and the corresponding power of the displacements can be determined analytically 
keeping trace of the singular terms in the massive Dirac equation. Only the exact 
constant in front of this singular term constitutes a problem, and it can be evaluated 
numerically. We shall address this case in detail elsewhere, and here only consider 
a few limiting cases. 

(i) Heavy and light quarks in the final state, e.g. b — > cdu 

It is possible to see that, in this case, the exponential terms, in the leading 
approximation, are associated with the light quark, and the presence of the heavy 
quark in the final state has no impact apart from changing kinematics. Indeed, 
we saw that the Green functions of the final state particles enter at large distances 
\x\ ~ p (it will be also illustrated in more detail in the next section). In this 
situation the interaction of the final heavy quark with the soft background field 
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reduces to the ordered exponential of i j A^)d^ along the heavy quark trajectory. 
We have already calculated it and found to be unity (the exponent to vanish) in the 
saddle point configuration. Therefore the final heavy quark can be taken as non- 
interacting. The propagator of the non-interacting final c quark oc e~ imc ' x ' (e~ mc \ x \ 
in the Euclidean time) will be multiplied by the light quark Green function, which 
develops a pole in the complex x 2 plane. This means that as far as the exponential 
terms are concerned, the effect of the final heavy quark is the replacement of wig in 
the Fourier transform by mq — m c , 

m Q — > A = m Q - m fin (53) 

(in the general case, with a few heavy particles in the final state, A = itlq — m &n > 
with the sum running over all heavy final state particles). Calculating the width, one 
thus has, at the saddle point, the exponent of the form e _2Ap rather than e~ 2mQP . It 
is worth emphasizing that this result holds as long as the mass of the final state quark 
involved exceeds 1/p, i.e. m fin p ^> 1, regardless of the actual hierarchy between mq 
and mfi n . The very same kinematic change refers, of course, also to the case when 
massive leptons (or other color-singlet particles) are present in the final state. 

Such a result might seem counter-intuitive to the reader who would compare 
it, say, with the free quark answer for the total width, where the final quark mass 
appears only as a quadratic correction oc m\ n /mQ if m^/mq -C 1, and the final 
heavy quark is typically fast. As a matter of fact, there is no mystery - the oc- 
currence of e~ 2Ap , instead of e^™^, is a manifestation of a remarkable property 
of the exponential terms discussed in Sect. 5. These terms are determined by the 
highest possible invariant mass of the light quark system. In the case at hand this is 
achieved when the final c quark is at rest, i.e. we are in the SV limit With this 
picture, it is then no surprise that the exponential terms are determined by e _2Ap . 

Since this result is rather unusual, we reiterate. In calculating the exponential 
terms in the transitions of the type b — > cud we find ourselves in the situation of 



the heavy quark symmetry |34], pql . It is an extended symmetry, however, since it 
applies to arbitrary color structure of the weak vertices, e.g. even when color flows, 
say, from Q to the light quark q rather than from Q to the heavy final quark c 
(the standard heavy quark symmetry works only for the color singlet be currents). 
Similar consideration applies even to the decays like b — > ccs. 

(ii) SV limit. 

In the SV limit, when the final heavy quark is slow, we can calculate the ex- 
ponential terms analytically in a wider class of processes, e.g. in the nonleptonic 
transition b —>■ cud. 

This case can be treated as follows. The final heavy quark Green function is 
given by the P exponent together with the "mechanical" phase factor (we write it 
here in the Euclidean space): 

1 i 7 '52 -m\x\ 

v ' ; 2 (27rm|x|) 3 / 2 
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The path exponent is again given by the U matrix presented in Eq. (^) []. At the 
saddle point when the instanton is on the line (x, y ) the final heavy quark path 
exponent equals to unity, and, therefore, the heavy quarks appear to be "sterile" - 
they do not feel the instanton field at all. Thus, in nonleptonic decays in the SV 
limit the appearance of exponential corrections is due - in the leading order - only to 
the light quark interaction with the instanton. (Note also, in the semileptonic decay 
in the SV limit the exponential terms are absent in the approximation accepted in 
the present paper. To see them we have to go beyond the leading approximation 
near the saddle point approximation. ) 

It is worth noting that, outside the exponential factor, in the pre-exponent, the 
masses enter in a more complicated way than is indicated in Eq. ([53]) . For example, 
the integration over d 3 x near the saddle point produces tuq 3 ^ 2 rather than the power 
of the energy release. 



6.3 Summary 

Now we are ready to incorporate the effect of both massive and massless quarks in 
the final state. 

1) Each light (massless) quark contains a second order pole in its propagator 
([43]) and brings in two powers of energy release in the numerator. 



2) Each heavy (static) quark propagator contains factor rri 



3/2 
fin 



3) The scale of energy release in the exponent is set up by the masses of the 
heavy quarks, 

rriQ - VJm nn 
Exploiting Eq. (|49|) and discussion above it we obtain 

T oc const • (m Q - m^^ 12 m Q 3/2 e ~ 2p ^-^ m ^ J]™! (55) 



where n\ is the number of massless quarks; iriq is the mass of the initial quark and 
m nn are the masses of the final heavy particles (both quarks and leptons). If there 
are no heavy particles in the final state we return to Eq. (49). 



7 Heavy quarks — a toy model 

Before proceeding to the actual calculation of the instanton contribution to heavy 
quark decays, we will first perform the same analysis in the simple toy model of 

6 Here, we can first deform the integration contour over xq into the complex plane and then use 
the l/mfi n expansion for the quark propagator. Another clarifying remark: for massive particles, 
the propagator entering at the saddle (effectively Minkowski) point is complex. In doing the saddle 
point calculations, one should take its values at the "bottom" of the cut; we shall dwell on this 
point elsewhere. 
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scalar chromodynamics. Without any loss of physical content, we can simplify our 
consideration by neglecting the spin degrees of freedom, as was suggested in Ref. 
f36fl . We also assume the SU{2) color group here. 
The weak Lagrangian, 

Lw = hQqcf) + h.c, (56) 

describes the decay of a heavy scalar quark, Q, into a light (massless) quark, q„ and 
a scalar "photon" , ; the coupling, h, has dimensions of mass. Both quarks are in 
the spinor representation of the color group. 

The basic strategy of our semiclassical calculation of the transition amplitude 
has been outlined in the previous section; here we work out details. We calculate 
the transition operator in Euclidean space using the semiclassical approximation, 
considering scalar quarks in the background of an instanton field. Upon returning 
to the Minkowski space, T will acquire an imaginary part related to the instanton 
correction to the total decay width. 

7.1 Inclusive width 

Consider the transition amplitude, 

T(r) = ^^( H Q\ f (r)\H Q ) = ^-(H Q \ J Ae "*™ iT{L w (x) L w (0)}\H Q ) 

(57) 

where r is an arbitrary external momentum. Denote the 4-velocity of the heavy 
hadron by (v = 0). Proceeding to the nonrelativistic fields, Q, and using Eq. (j38f), 
we have 

f(r) = i [ Q(x)S(x,0)Q(0)G <t> (x 2 ) e i{m ^~ r)x d 4 x d^z , (58) 

where is the propagator of the scalar photon, and S(x, y) is the propagator of the 
massless scalar quark in the external (instanton) field. We put y = 0. Addressing 
the total width, we will only consider r = and ro arbitrary (and complex); the 
more general case is relevant for differential distributions. For heavy quarks the 
transition amplitude depends only on the combination mg — tq : 

T(r) = T(m<g - r ) . 

Choosing an appropriate r , we select deep Euclidean kinematics and calculate the 
amplitude in the presence of an instanton. To this end we write 

m Q -r = ik , (x , z ) -> -i(x ,z ), (59) 

and 

T(ifeo) = J Q(x)S(x, 0)0(0) G+(x) e ikx d A xd 4 z, (60) 

where now everything is in Euclidean space; in what follows it is assumed that 
k ~ mg. 
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The propagator of the massless scalar particle in the instanton background is 



known exactly, and in the singular gauge takes the following form |33| 

where, to simplify the expression, we have assumed that the instanton is centered 
at z = and lies in a particular SU(2) color subgroup, 

A (t) - 2r l^ a ( X ~ Z )vP 2 f fi n\ 

A ^ X) ~ ( X - z )2(( x _ z )2 + p 2) ^ 

and 

r = (r, i) r+ = (f , -i) r+r^ = 5 a/3 + r/ Q/3c r c , (63) 

where f are the Pauli matrices acting in the color subgroup. 

The heavy quark field, Q(x), is the solution of the equation of motion in the 
instanton background 

iD Q(x) = -^—(iD) 2 Q(x) 
Iran 



which, in the leading order in 1/mn, yields Eq. (HSJ) with the matrices U given by 
Eq. (f46"|) . For the heavy scalar quarks one has 

' (H Q \QQ\H Q ) = -L . (64) 



2M Hq ^ 1 1 2m Q 

Now we take a closer look at the final state quark propagator, rewriting it using 
the Feynman parametrization: 

1 -V 1 x 



s(x, v ) = - f di -or 

7T J 



i(l - i){x - y) 2 + p 2 + P 



1 / p 2 r + (x — z)r(y — -z) 



1 + V(* - *) 2 V% - *) 2 (65) 



4ir 2 (x — y) 2 \ (x — z) 2 (y — zf 

where 

z = z-x£-(l-(i)y . 
Integrating in Eq. ( pD| ) over xo we only pick up the pole at 

2 P 2 + Z 2 

x — 



which yields the following exponential factor: 
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At k p ~ ttiqp ^> 1, the remaining integrations are nearly Gaussian, and run over 
narrow intervals, 



12 P (r 1 V . 1 ( . X Y . P 



x*~JL- , (67) 
m Q V 2/ m Q p V 2/ m Q 

Thus, one performs the remaining integrations by merely evaluating all pre expo- 
nential factors at the saddle point. In particular, this refers to the matrix U~ 1 (x), 
coming from the heavy quark propagation, which is now the path exponent from 
the point (x , ) to (0, ) . U^i^x) evaluated at the saddle point and is just the unit 
matrix, and the color part of the light quark Green function Eq. (^l|) . The heavy 



quark field, Q(0, x), according to Eq. (|67D, enters at distances x ~ J pj niq <C p and, 
therefore, the transition operator is finally proportional to Q(0) Q(0) ■ 
Collecting all remaining factors, one has 

f(Jfe ) = h 2 0(0) i ^" 4 ^^ / d£ d*xd A z e - fc oV(p 2 +- 2 )/«(i-0)+- 2 )l Q( ) , (68) 

2tv 2 p J 



where 



GJx 2 



4:7l 2 X 2 

is the free scalar propagator. Performing the Gaussian integrations (which yield the 
factor 47r 4 p 3 /cQ 4 ■ e~ 2fcoP ), and using Eq. fl6"4|) , we finally arrive at 



U2 -2k p 2i(niQ-r )p 

r( * ; = -«s^ = ~ h2 iwW^r ' (69) 



Let us note that the above equation shows the correct pre-exponential power of tuq. 
The propagator of the scalar particle in the instanton background has a first order 
pole, unlike the second order one of the spin-half quark, and matrix element fl64|) 
has an additional power of mQ 1 . This leads to a result with two powers of rriQp less 
then predicted in Eq. (0) . 

Now, we are interested in the decay width, which is given by the imaginary part 
of T(r) at the physical Minkowski point r = 0. Since the singularity of the amplitude 
we calculated is located far enough away, at r — rriQ, we perform a straightforward 
analytic continuation by merely setting r = and, thus, get 

In accordance with the general analysis of Sect. 3 the instanton contribution in the 
width decreases only as a power of tuq and oscillates. 
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7.2 Decay distribution 

Here we address the calculation of the instanton-induced corrections to the "photon" 
spectrum. The corrections blow up near the tree-level endpoint, where the semi- 
classical approach is not applicable. On the other hand, the total width is given by 
the integral over the whole spectrum and is calculable. The situation here is similar 
to the one we discussed in the Sect. 5 for the finite energy moment integrals in the 
e + e~ cross section, and we closely follow this analogy in our analysis. 

Let us denote the photon momentum by q, with q$ = E 1 . To find dT/dE^, one 
must consider the transition operator similar to the total width, but without the 
photon propagator : 

dT{r > q) = 2^r {HQldf(q ' r)lHQ) (71) 

Q 

with 

df(q,r) = i I Q(x)S(x,0)Q(0) e t(m o v - q - r)x d 4 x d 4 z . (72) 

The transition operator then clearly depends only on the sum of the four-momenta 
q + r, but since calculating the differential decay rate we need to keep q 2 = and 
qo = E y > 0, such a temporary proliferation formally allows us to use the momentum 
r for analytic continuation and is appropriate. In particular, we will keep q fixed, 
say, directed along the z axis, assume r = 0, and again use the variable tq to 
make calculations in the Euclidean domain. Kinematically, the amplitude dT(r, q) 
depends on two invariants, (r + q) 2 and (r + g)o, or any two combinations thereof. 

Technically, the calculation of the transition operator does not differ from the 
case of the total width except for the fact that now the Euclidean momentum k 

m Q - <?o ~ r o = ik , q + r = k 
has not only the zeroth, but also spacelike components. The saddle point calcu- 



lation goes exactly in the same way if one replaces k by \fk 2 = yk^ + k 2 , the 
main contribution comes from the singularity of the light quark propagator and the 
heavy quark matrix U still equals unity at the saddle point. For this reason, in the 
leading approximation the amplitude dT(r, q) appears to depend, in fact, only on 
one kinematic variable {rnqv — r — q) 2 = (p — q) 2 , and one has 

U2 2 -2pv / fe2 , 2 2 2ipy/(m Q v-q) 2 

dT(q,r) = n — — = n — — - . (73) 

m Q k 4 m Q (rriQV - g) 4 p 4 

In the last equation, we continued the result to the physical domain setting r = 0; 
to be far enough from the singularity and ensure the applicability of the calculations 
we must assume that (itiqV — q) 2 p 2 ^> 1. 

The differential decay rate for the massless is given by 

£ = ^o%o)ImrfT(g) = 
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h 2 o 2 sin (2p x /(m v - q) 2 ) 

= . » (mQ - 2,o) l ( ,^_; )4 ;' J (74) 

(remember, we fixed |g | = go)- Eq. (|73| ) (and fl74l ) for the imaginary part), in fact, are 
nothing but the soft instanton contribution to the massless (scalar) quark propagator 
(in the Fock-Schwinger gauge, or dressed with the path exponent to make it gauge- 
invariant in the general case). 

In the approximation of free quark decay, the photon spectrum is monochromatic, 
oc 5(E 1 — mq/2). Equation ( |7"1|) yields a decay spectrum below the end point, 
at E 1 < (mq/2). The result for the spectrum is, as expected, oscillating (sign- 
alternating). This does not lead to any physical problems, of course, because this 
contribution is to be considered on the background of "normal" OPE corrections 
(as well as the perturbative ones) which populate the spectrum below the two-body 
endpoint. The main OPE contribution near the end point is related to certain initial- 
state interactions of the heavy quark, and is interpreted as Fermi motion |32], |36], 37j. 
It produces a decay distribution which is unsuppressed, of order unity, in the interval 
\E 7 — i71q/2\ ~ Aqcd, and decreases fast only when mq/2 — E 1 becomes larger than 
a hadronic scale. At (mq/2) — E 1 ^> Aq C D the perturbative tail of the spectrum 
density takes over. In this domain our instanton calculations are already legitimate, 



since there ^{mqv — q) 2 ~ ^Jmqp if p ~ Aq C D- However, we emphasize that the 
corrections (|7"3|), (f?4|) are not related to the Fermi motion: the instanton effects in 
the latter appear in sub leading orders in mqp. 

Equation (|7|), far enough from the end point, represents, in a sense, a purely 
"exponential" effect. This is not the case anymore, however, if one attempts to 
integrate over E 7 and find the contribution to the total decay width. This fact 
is clearly revealed at the technical level when one compares the above calculation 
with the preceeding calculation of the total width: in the case of r to t one effectively 
uses the ^-independent photon propagator G^x 2 ) = — l/(16n 2 p 2 ), whereas Eq. ([T4]), 
according to Cutkowski's rules, corresponds to the propagator G^q) = —1/q 2 , i.e. 
G^x 2 ) = \/{Ati 2 x 2 ). They differ explicitly by the presence of the singularity at 
x = 0, which generates the normal OPE terms in r tot . Of course, a thoughtful 
instanton calculation of the integrated spectrum yields the same result as the direct 
calculation of the total width. Starting from the differential spectrum one has to 
resort to a special treatment of the OPE domain near the end point. 

A straightforward attempt to calculate the total width by integrating the spec- 
trum (|74]) seemingly faces a surprising problem: the correction grows fast toward 
the end point 0, [8| where the light quark is soft (though carrying large energy), and 
the integral seems to saturate at mq/2 — E 1 ~ l/(mqp 2 ), where the expansion in 
1 / Qp fails, and the overall result completely depends on an ad hoc cutoff procedure; 
the total contribution then would allegedly be governed by this soft scale, rather 
than by mq. This is similar to the apparent paradox in the finite energy moments 
of R in the e + e~ annihilation addressed in Sect. 5. The resolution of the paradox is 
essentially the same: any exponential contribution in the hard part of the spectrum 
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is automatically accompanied by the corresponding OPE-like "counterterms" in the 
end point region, which are superficially invisible. Although the exact spectrum in 
the latter cannot be calculated, certain integrals are defined unambiguously. As a 
result, the non-zero total width effect emerges only from the explicit (non analytic) 
constraint on the photon energy at E 1 = 0, i.e. at small q, and, thus, the instanton 
(exponential) contribution to the total width is indeed determined by the hard scale 
m Q- 

To reveal this conspiracy explicitly, we again exploit the fact that dT exponen- 
tially decreases in the Euclidean domain. In our concrete case, there is a technical 
simplification: dT(q) depends only on one kinematic variable k 2 = (ttiqV — q) 2 , the 
"invariant mass squared" of the light final state quark. One can write the dispersion 
relation over k 2 , 

9s 1 f 00 lmdT(K' 2 ) , ;9 . , 

dT{n 2 ) = - / , 2 2 dn' 2 . (75) 

7T JO K — K 

The fact that dT{n 2 ) falls off exponentially in the Euclidean domain, dT(n 2 ) ~ 
e -2V-K 2 p^ means that all moments of Im dT(n 2 ) vanish: 

/ K 2n lmdT(n 2 ) dn 2 = . (76) 
J o 

On the other hand, lm.dT{n 2 ) differs from dT /dqo by only a simple kinematic factor. 
For a massless 0, for example, 

= go %o) Im dT(g) , (77) 
dqo 2tt z 

rrin — k 2 



% = , Q = Qo 



Therefore, 



Mh ° 12 ■ , ''I 1 r frn 2 Q -KX'\, rri ,, 



dq g" — = / — y Im dTU 2 ) dn 2 . (78) 

dq 4:7T 2 m Q J m 2 Q -2m Q E \ 2m Q J 

We see that the instanton correction to the spectrum integrated from the very end 
point down to some E min is determined, not by the end point effects, but by the 



lowest energy included, and one has in this case 2pQ — > 2p s jrr?q — 2mQE m \ n . If one 
integrates over the whole spectrum, the correction is parametrically minimal, and 
reproduces the correction for the total width. Clearly, it is the most general feature 
of the "exponential" effects which does not depend on the particular decay. As was 
mentioned in Sect. 5, the constraints (|76|) can be written in the following compact 

form en 

lmdT{K 2 ) = / ds X {s) \5{s-K 2 )-e- s ^5(K 2 ) , (79) 
Jo 1 J 
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where x( s ) is some smooth function vanishing at s < 0. This expression automati- 
cally satisfies the relation 

rn 2 poo 

/ Im dT(n 2 ) P{k 2 ) dn 2 = - Im dTU 2 ) P(k 2 ) dn 2 (80) 

JO Jk 2 

for any appropriate analytic P(k 2 ). The saddle point instanton calculation carried 
out above determins the asymptotics of the function x( s ) at s 3> 1/p 2 : 



X (s) ^ n 2 h 2 ^- Sin(2 f P) . (81) 
m Q s 2 



In such a saddle point way we cannot calculate Im dT(n 2 ) at Kp < 1. According 
to Eq. ([F^), however, we know certain integrals over the small virtuality domain. 
For any particular "exact" xi s ) they are obtained as a series of 5 functions and 
derivatives of 6 functions located at the end point, which, when summed, give a 
certain more or less smooth function in the whole range, whose asymptotics are 
given by Eq. ( |8TD . 

For example, one can calculate the instanton contribution to the total width by 
integrating the spectrum. We use the general relation of the type (|76|) or (f78|) to 
write the total width as minus the integral of the right- hand side of Eq. ( |74|) (without 
step-functions) from q = — oo to 0. In the leading order in mgp, the integral is 
easily performed by parts, and is determined near the vicinity of the point go = 0, 
immediately yielding, exactly, the total width flTO]). 

It is not difficult to see how this works in the most general case. Suppose we 
study, for example, the "semileptonic decays" Q — > q + £ + v. If we can measure 
momenta of both I and u, pe, u , the differential distribution is merely given by 

d 2 T = <5>{ Pe ,p u )-lmdT{m Q v-{p e +p„)) . (82) 

Now we want to integrate over the momentum of neutrino and determine the in- 
stanton contribution to the charged lepton spectrum. We then keep pi fixed and 
integrate Im dT with the phase space, depending on p u . If the phase space fac- 



tor were an analytic function of p u , the constraints (|76|) or ([79]) would ensure the 
vanishing of the integral. 

However, the phase space contains a step-function at E u = 0, and only for this 
reason does one get a nonvanishing integrated width due to exponential terms. Then 
the relation (P0|) can be used to represent the integrated effect as the effect from the 
kinematic boundary and above (i.e., over the domain where neutrino "carries away" 
negative energy), which is a hard domain where the expansion is applicable []. Since 

7 The same reasoning equally applies, of course, to the case of the massive final particles as well. 
In this case it is convenient to phrase this consideration in the frame where vciqv — pi has only a 
timelike component. The phase space integral is not a polynomial anymore when m„ ^ 0, but is 
the step-function at E v = m v multiplying a fractional power function. Then one merely must put 
the corresponding phase space factor P in Eq. (|8^) under the sign of Im in the right hand side. 
This does not change the general reasoning presented here. 
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in this domain Im dT(n 2 ) is a rapidly oscillating function ~ sin(2fi;p), the integral 
is determined by its lower limit; the expansion in 1/Q is obtained by integrating 
by parts, and in the leading approximation given by the corresponding derivative 
of Im dT(K 2 ), from which we must keep only the derivatives of sin(2ftp). The 
exact coefficients combine to yield just the value of the neutrino propagator at the 
point which enters in the direct calculation of the lepton spectrum, i.e. G„(x) at 
x = 2p(m Q v - Pi)/\m Q v - p e \ . 

To illustrate the last assertion, we must remember that to get the total width in 
the leading approximation, we need to integrate the differential width near p v = 
keeping trace of only the non-analytic part due to the neutrino phase space, and 
the oscillating part, sin (2pJ {rrigv — pp — pv) 2 )., m ^ ne differential distribution. All 
other constituents of the amplitude can be approximated by their values at p v = 0. 

To see that this integration automatically yields the proper factor, we can use 
the following trick: compare this integral with the calculation of the neutrino Green 
function in the coordinate space by taking the Fourier transform of its momentum 
representation. Let us take the Fourier integral by closing the contour of integration 
over ko by its physical residue (we assume here that the zeroth coordinate coincides 
with the direction of the vector ttiqV — pe) : 

r d A k r d 3 k 

G ^> = / wri GAk) e " fa = / wm GAk) e ,fa • (83) 

The large-x asymptotics of the last integral is determined by the behavior of the 
phase space factor d 3 k/ko at small \k\ where it is non analytic. This factor, on the 
other hand, is exactly the same as in the double distribution if one identifies k in 
the above calculation with p u (both are nothing but d 4 k 8+{k 2 — m 2 u ) / (27r 4 )) . If one 
chooses the coordinate x in such a way as to have the same oscillating exponent 
e ikx on the right hand side of Eq. QS3|) as in dT ~ e 2 W Vq^-jv) 2 ; identifying 
k with p v , the values of the integrals will also coincide. The last condition just 
fixes the above stated value of the coordinate. This matching is rather obvious 
since the oscillating factor in dT came originally from evaluating the exponential 
e i(m Q v-p t -p v )x a ^ the saddle point x, which enters the neutrino Green function. 

Therefore, integrating over the neutrino momentum in the decay, we recovered 
our direct recipe of calculating the lepton spectrum, namely considering the problem 
as a two-body one, but using dT given by the product of the quark propagator in the 
instanton field, and the neutrino propagator at the saddle point x 2 = — 4p 2 . Since 
the integrated effect comes from the zero momentum of neutrino, it is governed by 
the same exponential (oscillating) factor determined by mq and p, \ mqv —pi\. We 
can repeat the very same consideration once more integrating now over the energy of 
the charged lepton; the leading term in the integral again comes from pi = 0, where 
the phase space of the lepton is non-analytic, the nearby integration reproduces 
Gp{— 2ip, 0), the argument of the neutrino Green function becomes the same, and 
we arrive at the total width obtained in the direct way in the preceeding subsection. 
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8 Instanton Contribution to Heavy Quark Decay 
Rates - Real QCD 

We now proceed to the case of actual b and c quarks; the modification required here 
is accounting for the quark spins. In principle, the analysis goes along the same 
lines as for the scalar quarks. However, since we deal with instantons - topologically- 
nontrivial configurations of the gauge field - the massless quarks acquire zero modes. 
They manifest the intervention of the infrared, long-distance effects in the presence 
of the instantons. This is an obvious defect of our simplified one-instanton ansatz. If 
we used the topologically trivial (but nonperturbative) configurations, the problems 
with the zero modes would be absent. 

Technically, this problem emerges already at the very first step: the Green func- 
tion of the massless quark is not defined in the field of one instanton since the Dirac 
operator has a zero mode. In order to perform estimates similar to the ones de- 
scribed in the previous sections, we need to regularize the massless quark Green 
functions in the infrared. We do it in the most naive way: introducing a small mass 
term m q . Now, the Green functions are well-defined but they have terms which 
behave like l/m q , which show up wherever the chirality-flip quark amplitudes oc- 
cur. However, when the weak interactions of the quarks are purely left-handed, the 
problem disappears since the zero modes do not contribute, and we can put m q = 
in the end. We are aware, of course, that this procedure is not fully self-consistent, 
but, hopefully, it works satisfactorily for our limited purposes - revealing the correct 
exponential dependence, as well as the power of mg, in the pre-exponent. It seems 
quite plausible that only the overall numerical coefficient will be modified in a more 
accurate analysis. 

Otherwise, the calculations go with minimal modifications. Let us outline the 
treatment of the semileptonic width. Once again, we write the transition operator 
(the integration over p will be restored at the very end) 

S~i2 I t / |2 _ 

T(r ) = Fl n Qql I Q{x)Y,S(x^)Y u Q^)L^(x) e^ v -^ x d^x d^z (84) 



where the weak polarization tensor for the lepton pair L„ v and the weak vertices are 
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Lp,(x) = - — (2x ll x u - x 5^) , r M = 7 M (l-7 B ) . (85) 

x° 

The heavy quark fields, Q, are nonrelativistic as in Eq. ([53), and are originally 
bispinors. However, solving the Dirac equation of motion in the limit mqp > 1 we 
get the fields in the form of Eq. (|45[) . The color matrix U has the same form as in 
Eq.©. 

The Green function of the light quark in the instanton background, S(x,y), is 



expanded in m q and has the following form |33 

1 



S{x,y) = — — Po{x,y) + G{x, y) + m q A(x, y) + 0(m 2 q ) 

Tfl q 
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where P is the projector on the zero modes (P does not contribute to Eq. (|34"D since 
it flips chirality.) Therefore, we can merely put m q = to arrive at 

-S{x,y) — - — = — - — G{x,y) — - — = 



2 v 2 2 

n\-l/2 ( 2 (T + C)(Tn)\ 



^(l + p 2 /C 2 )- 1/2 (l + p 2 M- 1/2 f^%(^0(rA)(r+ 7 )(rr ? )') 



47T 2 A 2 CV ' ^ V > r , , j \f? + 7l* 

where 

( = x — z , r] = y — z , A = x — y . 



S6) 



Now, the calculation differs from the case of the scalar quarks in minor technical 
details, namely, the different power of (x 2 + p 2 ) in the denominator, and the presence 
of the 7 matrices in the numerator. At the saddle point, {7 = 1, the heavy quark 
fields enter at the origin, the leptonic tensor must be evaluated at x* = (2ip, 0), and 
the final result is 



G 2 F\VQ q \ 2 



T ^ = -TTTTT^ -ge^-^ Q(0)t 7o PQ(0) (87) 

4:7T z p ti (i(m Q - r )) 

where P is the projector onto the instanton color SU(2) subgroup, P = (1, 1, 0)di ag - 
Since 

(H Q \Q l0 Q\H Q ) = 2M Hq 

we finally obtain 

T -^-IWi s ' m{2m " p) - (88) 

In terms of the free quark semileptonic width, 

G 2 F m%\V Qq \ 2 

1927T 3 ' 

we get 

2 967T 

r 'i = ~ r o o 7 sin ( 2m Qp) ■ ( 89 ) 

3 [jriQpf 

In a similar manner it is easy to find the expression for the differential distribu- 
tions in the semileptonic decays. We quote here the expression for the instanton- 
induced lepton spectrum (in its hard part, i.e. far enough from the end point) in 
the same approximation, 

m b dT l (b -> utv) _ 2 48vr 2 ( e\ cos (2m b p\/l - < 



2T dE t 3 (m b p) 5 V 2) (1 - efl 2 
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m c dT\c -> s(d)£u) 2 48tt 2 cos (2m ePx /r^ y 

2r dE t ' 3 (m cP )5 6 U Cj (1 - e )5/2 • ^ 

Here e = 2E^/mi > or e = 2E^/m c for the two decays, respectively. It is immediately 
seen explicitly, using the technique described in Sect. 7.2, that the integral over the 
spectrum reproduces the total semileptonic width ([88]). 

Now let us proceed to the nonleptonic decays with the massless quarks in the 
final state. Repeating the derivation above we get the non-leptonic width 

T 2 2 2567r 

r„ = r- r °i^^ sin(2mQ/7) (91) 

where c± are the standard color factors due to the hard gluons in the weak vertex 



For completeness we also give the instanton contribution to the inclusive radiative 
rate of the type b — > s + 7, 

r ' K l r °j^f sin{2m ^- (92) 

The above estimate refers to the yield of photons with all energies. In experiment one 
cuts off the low-energy photons, however. According to the previous discussion, the 
introduction of the lower cut off can change the estimate of the duality deviations. 
We will not submerge into further details regarding this effect here. 



9 Numerical estimates 

In this section, we present numerical estimates of the possible violations of the local 
duality in our model. The effects rapidly decrease with the energy release. They can 
be quite noticeable at intermediate energies, however. The inclusive decay width of 
the D meson, and is expected to be one of the prime suspects. Indeed, with the 
mass of the c quark only slightly over 1 GeV, one can expect sizable violations of 
duality. Another case of potential concern is the hadronic width of r. This width is 
also saturated at a similar mass scale. In these two cases, there exists at least some 
(quite incomplete, though) empiric information. It is natural to treat one of them 
as a reference point, in order to adjust the parameters of the model. Basically, we 
have only one such parameter, the overall normalization of the instanton density do, 
see Eq. (0). We will use the semileptonic D decays for this purpose. Then the 
second problem (r decays) can be used as a check that the model is qualitatively 
reasonable and does not lead to gross inaccuracies. As a matter of fact, this was 
already demonstrated in Sect. 5.1. Encouraged by this success we then take the risk 
to use the model for numerical estimates of the duality violating effects in various 
B decays. Although our model is admittedly imperfect, the numbers obtained can 
hopefully be viewed as valid order-of-magnitude estimates. 

8 We should note that here the subleading effects are probably significant: even within our simple 
saddle point calculation we formally had to discard, e.g. terms of the type 9!!/(2Qp) 5 ~ (4/(2Qp)) 5 . 
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9.i r sl (L>) 

This decay was numerically analyzed in the heavy quark expansion more than once 
||39||-||45||, with quite controversial conclusions. We will summarize here our current 
point of view f46|, deferring a brief discussion of the literature until Sect. 11. 



The parton semileptonic D decay width is given by 

T Q (D -> luX ad ) = ~ 1.03 -KT 13 GeV at m c = 1.35GeV, (93) 

where the strange quark mass is neglected. The comparison with the experimental 
value 

T exp (D -> luX) ~ 1.06 • 10~ 13 GeV (94) 

seems to be very good. However, there are corrections to the free quark estimate 
(p3|); both the perturbative and nonperturbative corrections calculated within l/m c 
expansion work together to noticeably decrease the theoretical prediction. 

(i) Perturbative corrections 

The one-loop perturbative correction to the width in the four-fermion decay is 
known since the mid-fifties [57]; for QCD one gets the factor r/ pert multiplying the 
theoretical formula for the width (i.e. T — » r ??pert)) 

2 / 9 25\ a s , s 

V P en = 1 - 3 (tt 2 - -) ^ . (95) 

This factor obviously decreases the theoretical prediction for the width, the question 
is how much. The answer for 7] peT t is not as obvious as it might seem, and depends 
on how m c is defined. 

Equation (f55|) implies that one uses the (one-loop) pole mass of the c quark in the 
inclusive rate. (As well-known, the notion of the pole mass is ill-defined theoretically 
J3L ^D], |5D| . It is safer to use the Euclidean mass, which also pumps away some of the 



a s corrections from the explicit correction factor ?7p ert . This decreases the mass and 
the coefficient in the correction simultaneously. The product m^.r) peTt is numerically 
stable, however, and for our limited purposes we can stick to Eq. (|95D and the 
one-loop pole mass.) The one-loop pole mass was numerically evaluated, say, in 
the charmonium sum rules, yielding [51J the number 1.35 GeV quoted above. This 
number is also in a good agreement with the heavy quark expansion for the difference 
m b - m c [0], 

m b - m c = M B -M D + A { -L - -L) + O (^) (96) 

where 

_ M BiD + Mb* D* 
m b,d - 1 • 
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Here, and in what follows, we use the notations 



(B\bja^b\B) 3 



and 



(B\b(iD) 2 b\B) 



" J = 2Mb ' (97) 



Substituting m^ olc ~ 4.83 ± .03 GeV J53L and a reasonable value of /i 2 (see below), 



in Eq. (|96"D, we again end up with m^ ole ~ 1.35 GeV. 

Controversial statements can be found in the literature concerning the value of 
/i 2 , (associated mainly with its different understanding) but the issue appears to be 
numerically unimportant for our purposes. 

Now, one has to establish the normalization point of a s in Eq. (p5|). Luke et 
al. suggested |54| exploiting the BLM prescription |55j for this purpose; it must be 



done in strict accord with the treatment of the mass. This leads to 

VpeTt « (1 - 0.25) (98) 

(for further details see Ref. jfj5|). This number turns out to be stable against the 
inclusion of the 0(a 2 ), and higher order corrections estimated in a certain approxi- 
mation 



(ii) Nonperturbative corrections 

Now let us examine the nonperturbative corrections. There are no corrections 
to the width that scale like 1 / m c , and the leading ones are given by the 1 / m 2 terms 

r,(^r , p »(i-M-i|). (99) 

While nfj is known, see above, /i 2 is not yet measured in experiment. We have to 
rely on theoretical arguments, which, unfortunately, are not completely settled yet. 
The original QCD sum rules estimate |57j yielded 

/4 = (0.5 ±0.1) GeV 2 . (100) 

We believe that the value /j, 2 ~ 0.5 GeV 2 is the most reasonable estimate available 
at present. It matches a general inequality, 131 



/4 > fi 2 G ~ 0.36 GeV , (101) 

and a more phenomenological estimate of Ref. 0]. Moreover, Eq. ( |100|) is 
marginally consistent with the first attempt of extracting fi 2 directly from data 
|6~l| , keeping in mind the theoretical uncertainties encountered there (other analyses 
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are in progress now). In any case, the effect of the kinetic operator on the semilep- 
tonic D width is modest, so that the impact of the uncertainties debated in the 
literature can lead to at most ~ 5% change in Y S \(D). 

Assembling all pieces together, numerically we get 

Tf h (D) ~ T (1 - 0.25 - 0.3 - 0.15) (102) 

where the corrections in the parentheses stand for the perturbative correction, the 
chromomagnetic and the kinetic energy terms, respectively. Thus, one is left with 



less than a half of the experimental width. According to Ref. |46] , the next order 
nonperturbative 0(1/ ml) effects apparently do not cure - and possibly deepen - 
the discrepancy. 

Pushing the numerical values of the parameters above, within uncertainties, to 
their extremes (but still, within acceptable limits), one can somewhat narrow the 
gap, but it is certainly impossible to eliminate it completely. Therefore, it is natural 
to conclude that the observed discrepancy in T S \(D), at the level of several dozen 
percent, is due to duality violations. 

We will make this bold assumption. The instanton contribution to r s i(.D), Eq. 
(H|), is then convoluted with the instanton density ( |l8l) to yield 

2 967r 

I%(D) = -T - d Q rg sin (2m c p ) . (103) 

Ignoring the sine on the right-hand side 0, and requiring this contribution to be 
0.5 T , we obtain 

d « 9 x 10~ 2 (104) 

(the values p = 1.15 GeV -1 and m c = 1.35 GeV are used). We will consistently 
exploit the above values of do and po in all numerical estimates in Sect. 9.2. 

If our approach is applied to the hadronic r width, deviation from duality comes 
out to be 



2 



r r (r -> hadrons) , An 

M^M"*^" 5 " 10 ' (105) 

i.e. quite reasonable. The fact that the numbers come out qualitatively reasonable 
in this case is also demonstrated by Fig. |l|. Let us note in passing that the 5% 
uncertainty in T(r — > hadrons) translates into ~ 30% uncertainty in the value of 
a s (m T ). 



9.2 Duality violation in b decays 

What is the magnitude of the anticipated effects in other situations where they are 
not yet determined experimentally, say, in beauty decays? Within our model the 

9 The value of sin (2m c po) is sensitive to how close the argument is to kir. This proximity is a 
very model-dependent feature, sensitive to small variations of parameters. Since we are aimed at 
conservative estimates all sine factors here and below will be consistently put equal to unity. 
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answer can be given. In what follows we will use expressions obtained in Sect. 8, 
convoluted with the instanton density QI8D , ignoring sines and cosines in numerical 
estimates. In this way we expect to get an upper bound on the duality violating 
contributions. This expectation is based on the following: (i) smearing with a 
more realistic finite-width instanton density will inevitably result in some extra 
(exponential) suppression, compared to the delta-function density; if the instanton 
density is distributed rather narrowly, the above suppression plays no role in D and 
r but will presumably show up at high energy releases characteristic to B decays; 
(ii) replacing sines and cosines by unity we increase the estimated value of the duality 
violations. 

Let us first address the simplest case, b — > ulv semileptonic width, where the 
energy release is the largest. To abstract as much as possible from the untrustworthy 
details of the instanton model, we can merely use the scaling behavior Eq. fl89|) . In 
other words, in all expressions below we keep only the pre-exponential factors. Then 

\V\b^ulv)\ 2 , 96vr , . , 

d ^ ^2xl(T J . (106) 



T (b^u£u) 3 (m b p o y 

A similar estimate for the radiative transition, b — > s + 7, based on Eq. (|92l), 
yields 

Next, we move on to processes with a heavy quark in the final state. Of partic- 
ular interest are duality violating effects in the Kobayashi-Maskawa allowed b — > c 
transitions. Consider first the semileptonic decays B — > X c Iv. As was discussed in 
Sect. 6.2, the instanton result for this process vanishes in the leading saddle-point 
approximation, while the subleading terms near the saddle point have not been cal- 
culated. To get an upper bound on the duality violations in the b c transition, in 
a rough approximation, we neglect all these subtleties and merely use our expression 
for the b — ► u replacing m& by nib — m c , 

\T\B^clv)\ 2 96tt 4 

J — -, t-t- ~ -«o^7 ^ ~ 3 x 10 . (108) 

Summarizing, the duality violating corrections are expected to be negligible in 
the semileptonic and radiative B decays, even in the b — >• c transitions. 

Let us proceed now to nonleptonic decays. Here, according to our model, the 
situation may somewhat change: deviations from duality jump up. Intuitively it is 
clear that the smallest effect is expected in the channel b — ► uud where the energy 
release is the largest. Specifically, 

\F l (b uud)\ 2 , „ 256tt „ . , 

~ - rf oc 2 _ AC ,_ : ^10" 3 , (109) 



r (6 — > uud) 3 45(m(,po) 
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where 

r (u -a\ it- G 2 F m b \ V ub\ 2 
T (b^uud) = 3r = ^ 

and 

c _ = c ; 2 = {a s (m b )/a s (M w )} 12 / 23 « 1.3. 

The effect further increases in the case of Kobayashi-Maskawa allowed nonlep- 
tonic transitions b — > cud. 

Considering the c quark as a static heavy quark which, according to section 6.2, 
does not interact with the instanton field at the saddle point, we obtain the following 
expression: 

\T l (b^cud)\ „, 7c 2 , + 3c 2 _ + 2c+c_ 

2dn — ± x 



Here 



T (b -> cud) 12 

16 7r 5//2 (m b — nic) 1 / 2 /m c \ 3 / 2 a 

'— «3xl0~ 3 . (110) 



5 oT W 



mi p 



r (6^curf) ~ 3r -0.5 = 0.5 



64tt 3 

with the factor 0.5 due to the kinematical suppression in the phase space associated 
with the c quark mass. The peculiar expression with the c± factors above emerges 
from the color matrices in the weak Lagrangian and in the quark Green functions 
after averaging over orientations of the SU(2) instanton over the color SU (3) group. 

However, one may worry that the c quark is not heavy enough since the parameter 
m c p ~ 2 is rather close to unity. Then, we may also try to consider another limiting 
case, and treat the c quark as massless, keeping m c only in the energy release. Then 
applying the technique from Sect. 6.1, we obtain 



|r I (6 -> cud) I 2 , „ 256tt 

— (i ~ 2 - d c 2 _ j- -r «8x 10" 3 , (111) 

r (o — > cud) 3 45po(m b — m c ' 4 



which is close, numerically, to the previous estimate. Note, the coefficients in front 
of c + and c_ are different since the heavy quark does not have color structure in the 
propagator while the massless quark does. 

The instanton contribution is enhanced by an order of magnitude due to the 
fact that the energy release is by a factor, (m b — m c )/m b ~ 0.7, smaller then in the 
channel with massless quarks, and due to the kinematical suppresion in the partonic 
width. 

Finally, let us discuss the duality violating contributions in the transition with 
two heavy quarks in the final state, b — > ccs, where they are believed to be the 
largest, for an obvious reason: the energy release is the smallest 0. This natural 
expectation does not contradict our model, although the enhancement we literally 



10 



Phenomenological analyses of the b — > ccs channel are presented, e.g. in recent works |62 
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get is rather modest, cf. Eqs. (|110|) and (|113|) . The SV approximation for the c 
quarks in the transition b — > ccs is justified; the simplest saddle point evaluation 
yields a non-vanishing effect due to the presence of the s quark in the final state. 
Performing the saddle point evaluation, one obtains 

r [b ->• ccs (d)) ~ d 372- — — sm (2p (m b - 2m c )J , 

7rm 6 (m;, — 2m c ) ,:1/ ' ! po 

(112) 

with 

IT/ 12 IT/ 12 „„5 

T (b-,ccs) ~ 3r -0.15= 0.15 - 1 g 4 [s 1 » , 

and 0.15 coming from phase space suppression due to two c quarks in the final state. 
Numerically, 



\T\b^ccs)\ 2c\ + c 2 _ 64 vr 



2 m 3 



7 d + -- — c - ^—r&6xl0- s . (113) 



T (b -> ccs) 3 (m b po) 5 {m b (m b - 2m c )) 3 / 2 



Concluding this section, let us reiterate our numerical findings. Using the saddle 
point approximation, and the instanton model to determine the nature of the sin- 
gularities in the quark propagators, in the complex (Euclidean) plane, we are able 
to derive some scaling relations for the duality violating effects induced by "soft 
background" fields. Keeping track of only the powers of relevant energies in the 
pre-exponent and neglecting the rest, we observe the following hierarchy: 

Ar 1 _ Ar 1 Ar 1 

— — (b — > ccs) ~ 2— — {b — > cud) ~ 5— — (b — > uud) 
To To To 

Ar 1 Ar 1 Ar 1 

~ 16 (b^clu) ~ 75 (b^sj) ~ 300 (b -»• ulv) . 

To r r 

Our numerical estimates are expected to be upper bounds within the particular 
mechanism of duality violations considered in the present paper. We hasten to 
add, though, that there exist physically distinct mechanisms, e.g. due to hard non- 
perturbative fluctuations, or those which may be somehow related to the spectator 
light quarks in the initial state, and so on. They deserve a special investigation. 



10 Drawbacks and Deficiencies of the Model 

Our instanton model of duality violations has obvious shortcomings. Although qual- 
itatively it correctly captures the essence of the phenomenon we want to model - 
transmission of a large external momentum through a soft background gluon field 
- the one-instanton ansatz itself is too rigid to be fully realistic. It has virtually 
one free parameter, the instanton size, and this is obviously not enough for perfectly 
successful phenomenology. The value of p we use is even somewhat smaller than that 
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usually accepted in the instanton liquid model |23[ . Correspondingly, the oscillation 
period comes out too large. For instance, in R(s), Eq. fl28|) , the oscillation period is 
almost four times as large as one observes experimentally. Figure ^| suggests that 
the typical oscillation length in i?(s) exp is ~ 0.6 GeV while from Fig.|I| we get 
y/s ~ 2.7 GeV. Even if we used the instanton liquid value, this would not narrow 
the discrepancy in any significant way. Moreover, the data seems to suggest that 
the oscillation length slowly varies as we move to higher energies. Our one-instanton 
ansatz is certainly incapable of reproducing this feature. The lesson we learn is that 
the soft background field has to be larger in scale and more sophisticated in shape. 

Another manifestation of the unwanted rigidity of the one-instanton background 
is the occurrence of zero modes for massless quarks. This phenomenon also leads to 
some inconsistencies in our treatment. For instance, the two-point function of the 
axial currents would not possess the necessary transversality properties. We essen- 
tially ignored this problem, keeping in mind the emphasis we place on qualitative 
aspects of duality violations. After all, our model is semi- quantitative, at best. 

The soft gluon fluctuations crucial in the duality violations are definitely not the 
ones constituting the dominant component of the vacuum. Indeed, if we used the 
instanton weight (p~8D, with d fitted to reproduce the duality violations in the D 
meson semileptonic decay, as the instanton density in the liquid model we would 
render this model disastrous. With our density in the instanton liquid model we 
would get the value of the gluon condensate ~ 25 times larger than it actually is. 

In addition to the above negative features, the model obviously misses other 
mechanisms which might also lead to violations of duality. The most noticeable is 
the absence of the impact of the initial light quarks in b decays. At the very least, 
one could suspect that they play a role in the formation of the soft gluon medium 
which, after the decay, has to transmit a large energy release. The influence of the 
initial light quarks would make the duality violating effects non-universal (they will 
be different, say, in mesons and baryons). At the moment we have no idea how 
to take into account this effect, nor we have any idea of how essential it might be 
numerically. 

On the positive side, we would like to stress again, that the model is general 
enough. One considers instantons only as a source of finite distance singularities 
in the quark Green functions, and for that purpose they may serve satisfactorily. 
Our procedure has very little to do with the full-scale instanton calculations of the 
type presented in Refs. p0| . In this sense, our calculations are much less 

vulnerable than the standard instanton exercises. The finite-distance singularities 
merely represent the mechanism of transmitting a large momentum through a large 
number of soft "lines" , with no hard lines involved (so that this mechanism does not 
appear in the practical OPE). Our point of view is pragmatic: experimental data 
clearly indicate duality violations, with an oscillating pattern, and so we reproduce 
this physical effect through fixed-size instantons. Eventually, comparison with data 
will lead to a better understanding of the relevant gluon field configurations and 
emergence of a model free from the drawbacks summarized above. 
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11 Comments on the Literature 



The present work intertwines many aspects of QCD - purely theoretical and phe- 
nomenological - in one junction; some of these aspects are quite controversial and 
cause heated debates. Therefore, it is in order to briefly review the literature in 
which the relevant issues were discussed previously 

We have already mentioned previous instanton calculations in -R e + e - and r decays 
|30| , |26| , |27| , |28| , |29| . Technically, they are very instructive and advanced. There are 
hardly any doubts, however, that the solitary instantons considered in these exercises 
do not represent, in the dynamical sense, typical relevant vacuum fluctuations. The 
fact that the corresponding estimates fell short of the experimentally observed effects 
is neither surprising nor frustrating. At the same time, the provocative suggestion of 



Ref. |T7| to use instantons for abstracting finite-distance singularities in the quark 



Green functions, was largely ignored. We try to develop this idea to its logical limits. 

Recently, the impact of the small-size instantons was analyzed in the spectral 
distributions of the inclusive heavy quark decays, within the formalism of HQET 
J?]. ||. The "part larger than the whole" paradox was first detected in these works: 
the instanton contribution to the spectra was found to be parametrically larger 
than the very same contribution to the decay rate. (The result was divergent at the 
boundaries of the phase space. This divergence is due to the fact that the isolated 
instanton density badly diverges at large p, and the instanton size is regulated by 
the external energy release, which vanishes at the end points.) The solution of this 
paradox was discussed at length above. The spectra near the boundaries of the 
phase space can not be calculated point-by-point in the present-day QCD. Still, 
the integrals over the spectra over a finite energy range touching the end-point are 
calculable. Integrating the instanton contribution, taking into account its peculiar 
analytical properties, automatically yields the effect which is determined by the 
far side of the smearing interval, rather than by the end point domain. The main 
subtlety lies in the process of separation of a "genuine" instanton contribution from 
the regular OPE. Our procedure automatically avoids double counting, an obvious 
virtue which is hard to achieve otherwise. 

The question of whether or not the semileptonic D decays are subject to no- 
ticeable duality violations is more controversial and is debated in the literature. 
Sometimes it is claimed that the OPE result ( PP| ) is compatible with experimental 
data, with no additional terms. The price paid is rather high, however: the mass 
of the c quark is then pushed up beyond 1.55 GeV (the value of /i^ is pushed down 
almost to zero). If the value of m c was that high, one would be in trouble in many 



other problems, e.g. the charmonium sum rules [51], the analysis of the bb threshold 



region ]5B| , and so on. Moreover, using a calculational scheme, which relies on large 
m c , leads to poor control of the perturbative series - a fact noted in Ref. [ 54| . We 



believe that the value of the product m b c rj peTt used above is realistic, which inevitably 
entails violations of duality in the D meson semileptonic decays in the ballpark of 
several units times 1CT 1 . Let us emphasize that our standpoint is testable experi- 
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mentally. One of possible tests is analyzing, say, the average lepton energy in the 
D meson semileptonic decays. This quantity is much less sensitive to the value of 
m c than the total width. Therefore, the task of detecting deviations from the OPE 
becomes much easier. An indirect proof may be provided by confirmation of duality 
violations in the r lepton rate in the ballpark of several units times 10 -2 . 

12 Conclusions and outlook 

At high energies the inclusive decay rates (e.g. r — ► u+ hadrons, or B — > X u + £v, 
or nonleptonic B decays) are represented by the sum of the transition probabilities 
into a very large number of possible final states. It looks like a miracle that these 
complicated sums, with various threshold factors, final state interactions and so on, 
reproduce a smooth quark (gluon) curve. This duality is explained by QCD. If the 
quark (gluon) cross section is calculated by virtue of the procedure known as the 
practical OPE, one expects that the inclusive hadronic cross section coincides with 
the quark (gluon) curve at large energy releases up to terms which are exponential in 
the Euclidean domain, and have a very peculiar oscillating pattern in the Minkowski 
domain, where they fall off relatively slowly, at least in the oscillation zone (Sect. 
5.2). Physically these exponential terms are associated with the transmission of large 
external momenta through the soft gluon medium. 

To model this mechanism mentioned above we suggest instanton-motivated es- 
timates. The instantons are used to abstract general features of the phenomenon 
and, to some extent, to gauge our expectations. We associate duality violations with 
the finite distance singularities in the quark Green functions due to soft background 
field configurations. The instanton-induced finite distance singularities produce a 
pattern of duality violations which closely resembles the indications (rather scarce, 
though) provided by current experimental data on e + e~ and r decays. The free 
parameters of the model are calibrated using these data. 

Let us examine, for instance, Fig. 0. which presents the differential hadronic 
mass distribution in r decays. From the first glance it is clear that significant (up to 
~ 20^-30%) violations of local duality are present in the whole accessible range. As 
a matter of fact, the "oscillating", duality- violating part of the effective V— A x V— A 
cross section, 5Ry-Ai is well approximated, in the region 0.8 GeV 2 - 3 GeV 2 , by the 
function 



The right-hand side is about —0.1 at s = 2.5 GeV 2 . If we take this function literally 
at all s, and reconstruct the corresponding correlator SU in the complex plane, we 
will find that it has no l/Q 2 expansion at all and, hence, would be omitted in any 
calculation based on the practical OPE. In the Euclidean domain, the corresponding 




(114) 
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-1 -0.5 0.5 1 1.5 2 2.5 3 

Figure 6: Experimental data for Rv-a( s )] the solid line at s > is a Bessel function 
(see text) drawn to guide the eye through the data points. The curve at the Eu- 
clidean Q 2 is the corresponding McDonald function. D(Q 2 ) = \{Q 2 dj 1 dQ 2 )Ii(Q 2 ) . 



H(Q 2 ) decreases exponentially, 

We plot R(s) corresponding to Eqs. ( |114| ) and |115| in Fig. ||. The value of 5U(Q 2 ) is 
extraordinarily small, ~ 3 x 10~ 4 , already at Q 2 = 1 GeV 2 . It reaches a "noticeable" 
30% level only at Q 2 = 0.1 GeV 2 , at the mass of the two pion threshold. And it is 
still (at least) as important as the usual perturbative corrections in the Minkowski 
domain, in the physical cross section, at as large of values of s as s ~ 2.5 GeV 2 ! With 
some reservations we can say that in the e + e~ annihilation and in the r decays we 
already have direct experimental evidence that such effects are significant. 

Less direct - but still quite convincing - arguments show that semileptonic charm 
decays also exhibit a similar phenomenon. There are good reasons to believe that 
such a situation is not exceptional. We argued that these effects represent a typical 
behavior of the strongly confining interacting theory in Minkowski space. 

Accepting our model, with all its drawbacks, for qualitative orientation we were 
able to achieve certain progress in relating various pieces of phenomenology to each 
other. First, we found that the observed 20 to 30% deviations in e + e~ annihilation 
and in the spectra of the r decays are consistent with significant corrections to the 
inclusive semileptonic D decay rate. At the same time, our estimates of duality 
violations rapidly decrease with increase of the energy scale, and produce seemingly 
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negligible effects in the inclusive decays of beauty. 

The strongest duality violations are expected to occur in the non-leptonic decays 
of B mesons, especially in those which contain two charmed quarks in the final state. 
In the transitions h — > ccs they are of order of 1%, while in the semileptonic and 
radiative decay rates deviations from duality fall off in magnitude to several units 
xlO" 4 . 

Although our estimates are universal - they do not distinguish, say, between B 
mesons and A& baryons - the model per se, taken seriously actually carries seeds 
of "spectator-dependency" . Indeed, the presence of extra light quarks in the initial 
state (baryons of the type A&) can help lift the chiral suppression of instantons, 
enhancing their weight compared to the meson case. In the case of A& the spectator 
quarks can naturally saturate the instanton zero modes for u and d quarks in the 
diagram incorporating both "Pauli Interference" and "Weak Scattering" -type pro- 
cesses. Then it is natural to expect stronger violations of duality. The argument 
is quite speculative, of course. This issue has not been investigated in detail. A 
dedicated analysis is clearly in order. 

The instanton model we suggest for estimating duality violations relies only on 
the most general features of instanton calculus, deliberately leaving aside concrete 
details. Instantons are taken merely as representatives of a strong coherent field 
configuration which have fixed size p > 1/Q, providing the quark Green functions 
with finite-distance singularities. In this situation we get a transparent picture 
of the corresponding duality violating phenomena. Technically, in the Minkowski 
kinematics the effect of the finite- x singularities can be viewed as an additional 
emission of a spurious particle, a "ghost", with an arbitrary mass n. The mass of 
the fictitious ghost has a smooth distribution, decreasing as some power of but 
oscillating in a more or less universal way. Equation (|114|) , with 7.5-\/i replaced by 
2pK, is an example. The peculiarity of the k distribution is a remarkable fact that the 
overall decay rate, with emission of the "ghost" , is always saturated at the maximal 
invariant mass of the ghost available in the process at hand. This qualitative picture 
can actually be converted into a kind of special diagrammatic technique for the 
ghost propagation. The instanton-motivated estimate of the strength of the "ghost 
coupling" is probably too crude. However, we think that such an approach, in a 
generalized form, may prove to be useful in describing violations of duality. 

Exploring duality violations is a notoriously difficult task. This field practically 
remains terra incognito, over the two decades since the advent of QCD. Our present 
attempt is only the first step. The model itself has an obvious potential for im- 
provement. Getting rid of the rigidity of the one- instanton ansatz, one may hope 
to achieve phenomenological success in describing fine structure of the duality vio- 
lating effects: the length of the oscillations and its modulations, and so on. (The 
pattern of experimental data suggests, perhaps, the presence of more than one scale 
of oscillations in the duality violating component. More accurate data that could 
provide more definite guidelines, are still absent.) 

To this end, it is necessary to consider as a background field a more generic 
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configuration, with more free parameters, say, an instanton molecule, or liquid type 
configurations. It is clear that the configurations relevant to the phenomenon under 
consideration have at least two scales built in, and one of them is significantly 
lower than 1 GeV. Such a project will require a lot of numerical work, however - an 
element we wanted to avoid at the first stage. 

Two other promising directions for explorations of duality violations are two- 
dimensional models and weakly coupled QCD in the Higgs phase. Both directions 
are much simpler than the actual QCD, and still the phenomenon is complicated 
enough so that the answer is not immediately clear. The first attempt of using 
the 't Hooft model for this purpose was made in Ref. [p4|l , which contains some 
initial observations. The potential of the model is clearly far from being exhausted. 
As for gauge theories with the spontaneous breaking of symmetry, calculation of 
the inclusive two-particle scattering near the sphaleron mass, revealing the typical 
pattern of the cross section, would be extremely instructive for QCD proper |25| . 



Let us note in this respect that lattice QCD, unfortunately, can add very lit- 
tle, if at all, to the solution of the problem of duality violations. The reason is 
quite obvious: all lattice simulations are done in the Euclidean domain, where all 
"exponential" terms are indeed exponentially small. Physically interesting and nu- 
merically important are these effects at large energies, deeply inside the Minkowski 
domain, and very far from the Euclidean domain where the lattice simulations are 
formulated. 
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^^Remember, these configurations definitely have very little to do with those determining the 
most essential features of the QCD vacuum. 



55 



References 



[1] K. Wilson, Phys. Rev. 179 (1969) 1499; 

K. Wilson and J. Kogut, Phys. Reports 12 (1974) 75. 

[2] V. Novikov, M. Shifman, A. Vainshtein and V. Zakharov, Nucl. Phys. B249 
(1985) 445. 

[3] N. Isgur and M.B. Wise, in B Decays, Ed. B Stone (World Scientific, Singa- 
pore, 1992); H. Georgi, in Perspectives in the Standard Model, Proceedings of 
the Theoretical Advanced Study Instituite, Boulder, Colorado, 1991, Ed. R.K. 
Ellis, C.T. Hill, and J. D. Lykken (World Scientific, Singapore, 1992); T. Man- 
nel, Heavy Quark Mass Expansions in QCD, in Proc. of the Workshop QCD 
- 20 Years Later, Eds. P. M. Zerwas and H.A. Kastrup (World Scientific, Sin- 
gapore, 1993), vol. 2, page 634; B. Grinstein, Ann. Rev. Nucl. Part. Sci. 42 
(1992) 101; M. Neubert, Phys. Reports 245 (1994) 259; M. Shifman, Lectures 
on Heavy Quarks in Quantum Chromodynamics, Preprint TPI-MINN-95-31-T 
[|hep-ph/951037711 . 

[4] E. Braaten, S. Narison and A. Pich, Nucl. Phys. B373 (1992) 581; F. Le 
Diberder and A. Pich, Phys. Lett, B286 (1992) 147; 
A. Pich, Nucl. Phys. Proc. Suppl. 39BC (1995) 326; 
G. Altarelli, P. Nason and G. Ridolfi, Z. Phys. C68 (1995) 257; 
P. Ball, M. Beneke, and V.M. Braun, Nucl. Phys. B452 (1995) 563; 
M. Neubert, Nucl. Phys. B463 (1996) 511; 



M. Shifman, Int. Journ. Mod. Phys. All (1996) 3195 ||hep-ph/9511469 . 



[5] M. Shifman, Theory of Preasymptotic Effects in Weak Inclusive Decays, in Proc. 
of the Workshop Continuous Advances in QCD, ed. A. Smilga, (World Scien- 
tific, Singapore, 1994), page 249 |[hep-ph / 9405246|1 ; Recent Progress in the Heavy 



Quark Theory, Talk at the V Int. Symp. on Particles, Strings and Cosmology 
- PASCOS - March 1995, Preprint TPI-MINN-95/15-T ||hep-ph/9505289|] . 



[6] H. Georgi and S. Osofsky, Nucl. Phys. B420 (1994) 94. 

[7] J. Chay and Soo-Jong Rey, Z. Phys. C68 (1995) 425; Z. Phys. C68 (1995) 431. 

[8] A. Falk and A. Kyatkin, Phys. Rev. D52 (1995) 5049. 

[9] B. Blok and T. Mannel, Preprint TTP-95-22 |hep-ph/9505288f . 



[10] V. Novikov, M. Shifman, A. Vainshtein and V. Zakharov, Nucl. Phys. B191 
(1981) 301. 

[11] E. Eichten and B. Hill, Phys. Lett. B234 (1990) 511; 
H. Georgi, Phys. Lett. B240 (1990) 447. 



56 



[12] J. Chay, H. Georgi, B. Grinstein, Phys. Lett. B247 (1990) 399. 

[13] I. Bigi, M. Shifman, N.G. Uraltsev, and A. Vainshtein, Phys. Rev. D52 (1995) 
196. 

[14] V. Novikov, M. Shifman, A. Vainshtein, and V. Zakharov, Fortsch. Phys. 32 
(1984) 585. 

[15] G. 't Hooft, in The Whys Of Subnuclear Physics, Erice 1977, ed. A. Zichichi 
(Plenum, New York, 1977), p. 943; 
B. Lautrup, Phys. Lett. 69B (1977) 109; 

G. Parisi, Phys. Lett. 76B (1978) 65; Nucl. Phys. B150 (1979) 163; 
A. Mueller, Nucl. Phys. B250 (1985) 327; 

for a recent review see A. H. Mueller, in Proc. Int. Conf. "QCD - 20 Years 
Later", Aachen 1992, eds. P. Zerwas and H. Kastrup, (World Scientific, Singa- 
pore, 1993), vol. 1, page 162. 

[16] M.S. Dubovikov and A.V. Smilga, Nucl. Phys. B185 (1981) 109 

[17] M.S. Dubovikov and A.V. Smilga, Yad. Fiz. 37 ( 1983) 984 [Sou. J. Nucl. Phys. 
37 (1983) 585]. 

[18] I. I. Balitsky and V. M. Braun, Phys. Rev. D47 (1993) 1879. 
[19] G. 't Hooft, Phys. Rev. D14 (1976) 3432. 

[20] C.G. Callan, R. Dashen, and D. Gross, Phys. Rev. D17 (1978) 2717. 

[21] M. Shifman, A. Vainshtein, and V. Zakharov, Nucl. Phys. B163 (1980) 46. 

[22] E. Shuryak, Nucl. Phys. B203 (1982) 93, 117; 

D. Dyakonov and V. Petrov, Nucl. Phys. B245 (1984) 259; ZhETF 89 (1985) 
361 [Sov. Phys. JETP 62 (1985) 204]; 

E. Shuryak, Nucl. Phys. B328 (1989) 102. 

[23] E. Shuryak, The QCD Vacuum, Hadrons and the Superdense Matter, (World 
Scientific, Singapore, 1988). 

[24] A. Zhitnitsky, Phys. Rev. D53 (1996) 5821. 

[25] M. Maggiore and M. Shifman, in Baryon Number Violation at the Electroweak 
Scale, Eds. L. Krauss and S.J. Rey, (World Scientific, Singapore, 1992), page 
153. 

[26] LI. Balitsky, Phys. Lett. B273 (1991) 282. 

[27] LI. Balitsky, M. Beneke, and V.M. Braun Phys. Lett. B318 (1993) 371. 



57 



[28] P. Nason and M. Porrati, Nucl. Phys. B421 (1994) 518. 

[29] S.V. Faleev and P.G. Silvestrov, Preprint Budker INP 95-73 [[hep-ph/9509231 



[30] R. Carlitz and C. Lee, Phys. Rev. D17 (1978) 3238; 
N. Andrei and D. Gross, Phys. Rev. D18, (1978) 469; 
C. Baulieu et al. Phys. Lett. B77 (1978) 290; 
T. Appelquist and R. Shankar, Phys. Rev. D18 (1978) 2952. 

[31] T. Coan et al, (CLEO Collaboration), Phys. Lett. B365 (1995) 580. 

[32] I. Bigi, M. Shifman, N. Uraltsev and A. Vainshtein, Phys. Rev. Lett. 71 (1993) 
496; 

I. Bigi, B. Blok, M. Shifman, N. Uraltsev and A. Vainshtein, The Fermilab 
Meeting Proc. of the 1992 DPF meeting of APS, Eds. C. H. Albright , P. 
Kasper, R. Raja and J. Yoh, (World Scientific, Singapore, 1993), vol. 1, page 
610. 



[33 
[34 

[35 

[36 

[37 

[38; 



[39 
[40 
[41 

[42; 
[43; 
[44; 
[45; 



L.S. Brown et al, Phys. Rev. D17 (1978) 1583. 

M. Voloshin and M. Shifman, Yad. Fiz. 47 (1988) 801 [Sov. J. Nucl. Phys. 47 
(1988) 511]. 

N. Isgur and M. Wise, Phys. Lett. B232 (1989) 113; Phys. Lett. B237 (1990) 
527. 

I. Bigi, M. Shifman, N.G. Uraltsev and A. Vainshtein, Int. Journ. Mod. Phys., 
A9 (1994) 2467. 

M. Neubert, Phys. Rev. D49 (1994) 4623. 

M.K. Gaillard and B.W. Lee, Phys. Rev. Lett. 33 (1974) 108; 
G. Altarelli and L. Maiani, Phys. Rev. B52 (1974) 351; 

for definitions see e.g. Fayyazuddin and Riazuddin, A Modern Introduction to 
Particle Physics, (World Scientific, Singapore, 1992), Sect. 15.3.4. 

M. Luke and M. Savage, Phys. Lett. B321 (1994) 88. 

I. Bigi and N. Uraltsev, Z. Phys. C62 (1994) 623. 

Z. Ligetti and Y. Nir, Phys. Rev. D49 (1994) 4331. 

V. Chernyak, Nucl. Phys. B457 (1996) 96. 

A. F. Falk, M. Luke, and M. J. Savage, Phys. Rev. D53 (1996) 6316. 
M. Shifman and N.G. Uraltsev, Int. Journ. Mod. Phys. A10 (1995) 4705. 
N.G. Uraltsev, Int. Journ. Mod. Phys. All (1996) 515. 



58 



[46] B. Blok, R.D. Dikeman, and M. Shifman, Phys. Rev. D51 (1995) 6167. 

[47] R.E. Behrends, R.J. Finkelstein, and A. Sirlin, Phys Rev. 101, (1956) 866; 
S. Berman, ibid 112, (1958) 267 . 

[48] I. Bigi, M. Shifman, N.G. Uraltsev and A. Vainshtein, Phys. Rev. D50 (1994) 
2234. 

[49] I. Bigi and N.G. Uraltsev, Phys. Lett. B 321 (1994) 412. 

[50] M. Beneke and V. Braun, Nucl. Phys. B426 (1994) 301. 

[51] V. Novikov et al, Phys. Reports 41 (1978) 1. 

[52] A. Falk and M. Neubert, Phys. Rev. D47 (1993) 2965; 2982. 

[53] M.B. Voloshin, Int. Journ. Mod. Phys. A10 (1995) 2865. 

[54] M. Luke, M. Savage, and M. Wise, Phys. Lett. B345 (1995) 301. 

[55] S. Brodsky, G. Lepage and P. Mackenzie, Phys. Rev. D28 (1983) 228. 

[56] P. Ball, M. Beneke, and V. Braun, Phys. Rev. D52 (1995) 3929. 

[57] P. Ball and V. Braun, Phys. Rev. D49 (1994) 2472; for an update see E. 
Bagan, P. Ball, V. Braun and P. Gosdzinsky, Phys. Lett. B342 (1995) 362; 
earlier estimates can be found in M. Neubert, Phys. Rev. D46 (1992) 1076; E. 
Bagan, P. Ball, V. Braun, and H. Dosch, Phys. Lett. B278 (1992) 457. 

[58] M. Voloshin, Surv. High En. Phys. 8 (1995) 27. 

[59] M. Shifman, N.G. Uraltsev and A. Vainshtein, Phys. Rev. D51 (1995) 2217. 

[60] I. Bigi, A. Grozin, M. Shifman, N.G. Uraltsev, and A. Vainshtein, Phys. Lett., 
B339 (1994) 160. 

[61] M. Gremm, A. Kapustin, Z. Ligeti, and M. B. Wise, Preprint CALT-68-2043 
||hep-ph/96033I41 1; 

V. Chernyak, Preprint Budker INP-96-23 ||hep-ph/9604376| . 

[62] A. Falk, M. Wise and I. Dunietz, Phys. Rev. D51 (1995) 1183; 

G. Buchalla, I. Dunietz, and H. Yamamoto, Phys. Lett. B364 (1995) 188. 



59 



